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Abstract. We study an integrable modification of the focusing nonlinear Schrodinger equation from the 
point of view of semiclassical asymptotics. In particular, (i) we establish several important consequences of 
the mixed-type limiting quasilinear system including the existence of maps that embed the limiting forms 
of both the focusing and defocusing nonlinear Schrodinger equations into the framework of a single limiting 
system for the modified equation, (ii) we obtain bounds for the location of discrete spectrum for the associated 
spectral problem that are particularly suited to the semiclassical limit and that generalize known results for 
the spectrum of the nonself adjoint Zakharov-Shabat spectral problem, and (iii) we present a multipaxameter 
family of initial data for which we solve the associated spectral problem in terms of special functions for all 
values of the semiclassical scaling parameter. We view our results as part of a broader project to analyze the 
semiclassical limit of the modified nonlinear Schrodinger equation via the noncommutativc steepest descent 
procedure of Deift and Zhou, and we also present a self-contained development of a Riemann-Hilbert problem 
of inverse scattering that differs from those given in the literature and that is well-adapted to semiclassical 
asymptotics. 



1. Introduction 

This paper will be concerned with the modified nonlinear Schrodinger (MNLS) equation for a complex- 
valued function <f>{x,t): 

The correct problem to pose for this equation when x G M is the Cauchy problem, where initial data 4>{x, 0) 
are given at t = 0. Here a and e are real parameters which we assume throughout the paper satisfy the 
inequalities given in ([1]). 

The focusing nonlinear Schrodinger (NLS) equation is a special case of ([T]) when a = 0. The focusing NLS 
equation is of course very well known, arising naturally in the modeling of any nearly monochromatic, weakly 
nonlinear, dispersive wave propagation. In particular, it is in many circumstances an adequate model to 
describe the evolution of the electric field envelope of picosecond pulses in monomode optical fibers [51] . The 
so-called "semiclassical" or zero-dispersion limit of the focusing NLS Cauchy problem consists of taking the 
initial condition in the form 0) ~ 74(x)e*'^^^-'/'^ for fixed real functions A{-) > and S{ ) and considering 
the behavior of the corresponding solution 0(a;, t) of the focusing NLS equation in the limit e J, 0. This 
asymptotic problem is important in the design of increasingly prevalent dispersion-shifted fibers as well as 
being of intrinsic mathematical interest. 

The utility of the focusing NLS equation as a mathematical model stems from the fact that this equation 
is integrable, and in principal the Cauchy problem may be solved using the scattering/inverse scattering 
transform [41] that is precisely adapted to this equation. However, in applications if higher-order effects 
become important then it is necessary to account for them in the model. To consider sub-picosecond 
electromagnetic pulses in monomode optical fibers, one should consider the more general evolution equation, 
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dt ' 2 dx^ ' ^ dx 
where a, a', and a" are real constants and where the three additional terms on the right hand side account 
for nonlinear dispersion, intrapulse Raman scattering, and higher-order linear dispersion respectively [1] . In 
general, ([2]) is not integrable; however it turns out that there are special cases of ([2]) besides a = a' = a" = 
that are integrable, although by different machinery than applies to the focusing NLS equation. In fact, the 
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MNLS equation ^ is such an integrable system [ID] (related to the so-called "derivative" NLS equation). As 
the MNLS equation takes into account at least the higher-order effect of nonlinear dispersion, understanding 
the predictions of this equation can serve as a mathematical "stepping stone" in the understanding of the 
more general problem ^ as well as any weakly nonlinear, nearly monochromatic, dispersive wave propagation 
problem when nonlinear dispersion is not negligible. 

One of the well-known features of the focusing NLS equation is the so-called modulational instability. 
This instability may be manifested in several different ways, and so the term "modulational instability" can 
mean different things in different contexts. However, the simplest context of modulational instability is the 
linearized perturbation theory of small disturbances of exact "plane-wave" solutions of the focusing NLS 
equation. Indeed, the NLS equation 
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(the focusing case is k = 
relation 
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: 1 and the defocusing case is k = —1) has exact solutions (plane waves) of the form 
for real constants A > k and lo necessarily related by the nonlinear dispersion 
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Consider a small perturbation of this plane-wave solution by substituting 
is a new unknown. In this way, one obtains an (equivalent) equation for the 
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in-phase" perturbation p: 
-:-^ + kA\p + p* + 2b|2 + p2 + |p|2p) ^ , (5) 



where the asterisk denotes complex conjugation. Assuming p to be small, this equation may be linearized 
and written as a system for a = Rc{p} and h ~ Im{p}: 
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This is a linear constant-coefficient system of partial differential equations that may be solved using Fourier 
techniques. We look for elementary solutions of the form 
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where A is a "relative wavenumber" (to k) , which reduces 
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to a system of ordinary differential equations: 
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There is a basis of solutions proportional to eigenvectors of the coefficient matrix on the right-hand side, and 
the time dependence enters through factors of the form e"^*/^ where a is the corresponding eigenvalue of the 
coefficient matrix. The eigenvalues are given by 

(7 = -ifcA±^V4«;A2- A2. (9) 

If we are in the defocusing case (k = —1) then <t is always purely imaginary, and all Fourier-type solutions 
of ([6]) arc oscillatory in t. This behavior indicates the absence of modulational instability in the defocusing 
case. However, if wc are in the focusing case (k — 1) then a has a nonzero real part if A2 < AA^ and therefore 
there is a band of unstable relative wavenumbers about A = (and hence a sideband of the unperturbed 
wavenumber k). This is the modulational instability in its simplest form. These calculations show that 
in the focusing case, the only perturbations that have a chance of avoiding growth correspond to waves of 
length 0(e), and worse yet, the exponential growth rates for the unstable modes scale as In this sense, 
the modulational instability is enhanced when e is small. 

The modulational instability may also be understood from a nonlinear perspective, via the so-called 
modulation equations. Without loss of generality, we write the solution (j}{x,t) of ([3|) in the form (f>{x,t) ~ 
A(2;,t)e'''5(^'*)/'^ where A{x,t) > is a real amplitude and S{x, t) is a real phase, and then ([3]) implies 
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corresponding (real) nonlinear equations governing these two functions. If we think of A{x, t) and S{x, t) 
as functions independent of e (this assumption is of course not exactly consistent with ([3])) then locally 
in the neighborhood of fixed x and t, (j){x,t) resembles a plane wave solution with amplitude A = A(x^t\ 
wavenumber k = dS/dx{x, t), and frequency lu — —dS/dt{x, t). This observation is the grounds for expecting 
a connection between the preceding perturbative calculation and what will now follow. To proceed, it is more 
natural to introduce u{x,t) := dS/dx and p{x,t) := and obtain equations for these unknowns equivalent 
to ©: 



This time, the approximation we make (in lieu of the formal linearization step in the preceding discussion) 
is to drop the formally small terms in e. Thus, one arrives at the modulation equations associated to ([3]), 
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A system of quasilinear partial differential equations of the form (|lip is called hyperbolic if the coefficient 
matrix of the x-derivatives has distinct real eigenvalues and elliptic if the eigenvalues are complex. Cauchy 
problems for hyperbolic systems can be solved by the method of characteristics, one implication of which is 
continuous dependence of the solution on initial data for small time. The hyperbolic Cauchy problem is said 
to be well-posed. On the other hand, the Cauchy problem is ill-posed for elliptic systems. Worse yet, the 
Cauchy problem for an elliptic system can only be solved at all (by the method of Cauchy-Kovalevskaya) 
if the initial data p(a;, 0) and u{x,0) are analytic functions of the variable x\ this should be regarded as a 
overly restrictive condition for a physically relevant mathematical model. 

It is easy to check that for u € M and p > 0, the system ([TT|) is elliptic for k = 1 and hyperbolic for 
K = — 1 . This calculation therefore reveals the same dichotomy between the focusing and defocusing versions 
of the NLS equation as did linearization about a plane wave. The linearization calculation may be viewed as 
a "local in a;" version of the deduction of type (hyperbolic or elliptic) of modulation equations. In general, 
hyperbolicity of plane- wave modulation equations is equivalent to the absence of linear instabilities of plane 
waves of sideband type, although it is possible for a system to have hyperbolic modulation equations while 
admitting instabilities of relatively short waves |33| . Interestingly, analyticity of initial data as required 
by elliptic modulation equations plays a fundamental role in the spectral and inverse-spectral analysis of 
the focusing NLS equation in the semiclassical limit [57j, even though the latter analysis is not based on 
modulation equations at all. (For an approach based directly on the Cauchy-Kovalevskaya series solution of 
the analytic data Cauchy problem for the elliptic modulation equations, see Gerard [12] ■) 

One might notice that in the linearization calculation the conclusion of stability or instability was inde- 
pendent of the pair (fc. A) that characterizes the underlying plane wave. Indeed, in the focusing case, there 
is a band of unstable relative wavenumbers no matter what values k and A take. Similarly, the modulation 
equations (|lip are elliptic or hyperbolic independently of u S M and p > 0. In principle, it could have been 
otherwise; some waves may be stable to all sideband perturbations while others are not. For such a system, 
the modulation equations may be hyperbolic for some values of u and p and elliptic for others, and the 
system is said to admit a change of type. It will turn out that the system of modulation equations for the 
MNLS equation admits a change of type. 

To summarize, the focusing NLS equation exhibits modulational instability, while the defocusing NLS 
equation is modulationally stable. Now, the MNLS equation ^ appears to be a perturbation of the focusing 
NLS equation (through the limit a ^ 0) . One might expect on these grounds that the MNLS problem should 
also experience modulational instability. On the other hand, consider making the following substitution in 




(10) 



0(:r,O=e^(^«+^'-/2)/-V^(^,r), 
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where = x — ct and t = t, and c G R is an arbitrary parameter. Then, by direct calculation, ip satisfies 

+ y + (1 - V + ^ae-^^m'i') = . (13) 

If a = 0, we recover the well-known result that the focusing NLS equation is invariant under the one- 
parameter group of Galilean transformations p2|) . More generally, the MNLS equation is certainly not 
invariant under this group. In particular, one may note that by choosing c so that ac > 1. the sign of 
the term becomes negativcf^, and thus the equation governing 'il}{^,T) takes the form of a perturbed 

defocusing NLS equation rather then a perturbation of the focusing NLS equation. 

Given that when a = one expects drastically different behavior in the semiclassical limit depending on 
whether the sign of the cubic term is positive (strong modulational instability) or negative (modulational 
stability); the above calculation suggests two things: (i) the term in the MNLS equation ([1]) proportional to 
a certainly does not constitute a small perturbation uniformly in the phase space of fields 4> (that is, taking 
a 7^ amounts to a singular perturbation of the NLS equation), and (ii) that part of the phase space of 
fields may evolve under the MNLS equation ([T]) in such a way as to avoid all catastrophic modulational 
instability associated with the undifferentiated cubic term If this latter statement can be placed on 

rigorous footing, there are immediate implications for the design of optical fiber telecommunication links (for 
example). Indeed, one of the difficulties arising in long-distance fiber-optic systems is the so-called Gordon- 
Haus effect, in which pulses experience a random "jitter" in their arrival times which leads to detection 
errors. This effect has been shown to be a consequence of the modulational instability associated with the 
focusing nonlinearity (and anomalous dispersion) in the focusing NLS model. The possibility of avoiding 
Gordon-Haus jitter in the context of the MNLS model has already been explored [18| at the level of single 
soliton solutions. One of our goals is to understand this change in stability for the MNLS equation at the level 
of more general initial data, and in particular initial data that "contain" many solitons as is characteristic 
of the semiclassical scaling. 

Due to its integrable structure, solving the Cauchy problem for the MNLS equation (jT]) requires three 
steps: 

Step 1. Calculate certain "spectral data" associated with the given initial condition for the problem, 
viewed here as a coefficient in a linear system of differential equations with a spectral parameter. 
This is a "forward transform" of the initial data. 

Step 2. Determine how the spectral data change as the field 4) evolves under the MNLS equation, 
and obtain the spectral data corresponding to the evolved field. 

Step 3. Find the field (t>{x, t) from the time-evolved spectral data, that is, invert the mapping from 
step 1. This is an "inverse transform" step that in many cases can be written as a matrix- valued 
Riemann-Hilbert problem from analytic function theory. 
The miracle of integrability is that for an appropriate linear system of differential equations with spectral 
parameter, the inverse of the spectral map of step 1 that is required in step 3 actually exists, and as 
importantly the evolution of the spectral data in step 2 is completely explicit. For semiclassical problems, 
the small parameter e enters significantly in the first and third steps. Therefore, to understand the stability 
transition in the MNLS equation and its role in semiclassical asymptotics, we must first understand the 
forward transform (step 1) in the semiclassical regime. Once the the semiclassical asymptotics of the forward 
transform have been understood, it will remain to utilize the growing library of tools of asymptotic analysis 
for Riemann-Hilbert problems to complete the solution of the Cauchy problem for the MNLS equation in 
the limit e i 0. In this paper, we will consider the spectral problem associated with the MNLS equation (see 
(|176p ) paying particular attention to results that are meaningful in the limit e | 0. 

This paper is organized as follows. In Section[2]wc follow the procedure applicable to the NLS equation to 
obtain the modulation equations associated to the MNLS equation. Through the analysis of these modulation 
equations we will obtain a condition for modulational stability. Throughout the remainder of the paper we 
note the importance of this condition and the role it plays in the spectral problem. In Section [3] we establish 
some bounds on the discrete eigenvalues of the spectral problem (|176p associated to the MNLS equation. 



If c is chosen so that ac = 
derivative NLS equation for '0. 



1 exactly, then under the Galilean boost the MNLS equation for <p reduces to the so-called 
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In Section 3] we use the theory of special functions to calculate the spectral data explicitly for a certain 
multiparameter family of initial conditions; this calculation is valid for all e > and provides an explicit 
avenue to the scmiclassical analysis of the corresponding Cauchy problem. The appendix contains the general 
spectral and inverse-spectral theory associated with the MNLS equation (that is, all theoretical details of 
steps 1-3 above) and in particular includes a derivation of the associated Riemann-Hilbert problem that we 
will return to in a later publication. The theory outlined in the appendix frames all of our analysis of the 
MNLS equation. 



Notation. Throughout the paper, we use boldface capital letters to refer to square matrices, with the 
exception of the identity matrix which we write as I, the Pauli matrices 
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and certain diagonal matrices denoted 6°°^^ for a G C and defined by 
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(14) 



(15) 



Lowercase boldface letters refer to column vectors. We use the "dagger" notation (superscript f ) to indicate 
the transpose and elementwise complex conjugate of a matrix or vector, and we use asterisk notation (su- 
perscript *) for complex conjugation. We sometimes use prime notation for derivatives of functions of just 
one variable. 



2. Modulation Equations 

In this section we derive the modulation equations associated to the MNLS equation and find a condition 
on the local amplitude and wavenumber equivalent to the hyperbolicity of these equations. By analogy with 
the connection between hyperbolicity of modulation equations and stability of plane waves described in the 
context of the NLS equation in Section [U this calculation gives a condition under which the MNLS equation 
is (locally) modulationally stable. The fact that one can have local stability without global stability is related 
to the fact the MNLS modulation equations admit change of type. 



2.1. Derivation of the modulation equations. Consider the MNLS equation ([T]), and for some real- 
valued functions A{x,t) and S{x,t) write the solution (f>(x,t) as: 



(x,i) =^(a;,t)e'^^'''*^/' 



(16) 



Substitution into ([T]), canceling the phase factor e^^^^, taking the imaginary part, multiplying by A and 
dividing by s gives: 
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Proceeding similarly but taking the real part instead, one arrives at: 
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Setting p{x,t) :~ A(x,t)'^ and u{x,t) :~ dS{x,t)/dx, the system of modulation equations is obtained by 
dropping the formally small terms in e: 
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2.2. Modulational stability criterion. The eigenvalues ( of the coefheient matrix of the x-derivatives in 
p9)l are the roots of the equation: 

- {Aap + 2u)C + 'Saup + ia^p^ + it^ + p = 0. (20) 

The discriminant of this equation is 4p(Q;^p + au — 1). Since p > 0, the condition 

a^p + au - 1 > 0, (21) 

is equivalent to having real distinct eigenvalues, which makes (|19p a hyperbolic system. Note that for a given 
value of the parameter a e M, this condition involves both dependent variables p (square of local plane- wave 
amplitude) and u (local wavenumber) and this dependence is characteristic of a quasilinear system that can 
change type. The inequality (|2ip will appear at every stage of our analysis of the spectral problem (see 
(|176p ) associated with the MNLS equation. 

2.3. Implications of the stability criterion on Riemann-Hilbert analysis of the inverse spectral 
problem. While this paper is concerned mostly with the semiclassical limit of the spectral transform as- 
sociated with the MNLS equation and a subsequent paper will analyze the corresponding inverse-spectral 
transform in this limit, it is worthwhile pointing out how we might expect the stability criterion (|2ip to 
influence the inverse-spectral part of the solution of the MNLS Cauchy problem for rapidly decaying initial 
data of the form cj){x, 0) = ^(a;)e*■^(^)/^ a; € R. 

How do modulated plane waves emerge from the asymptotic analysis of the inverse-spectral problem? 
As shown in the appendix, the inverse-spectral problem for the MNLS Cauchy problem can be cast as a 
Riemann-Hilbert problem in which one seeks a 2 x 2 matrix-valued unknown depending principally upon 
a complex variable k and parametrically upon x, i, a, and e. The unknown matrix is characterized by a 
normalization condition as k — > 0, prescribed discontinuities along the real and imaginary fc-axes, simple 
poles with prescribed residues at certain other points in the complex fc-plane, and analyticity in k for all other 
/c G C. Such a Riemann-Hilbert problem involving a parameter e in a singular way may be analyzed by a 
two-step process. First, one introduces an explicit transformation of the matrix-valued unknown designed to 
remove the poles at the expense of introducing new jump discontinuities. Thus one arrives at an equivalent 
Riemann-Hilbert problem for a piecewise-analytic matrix-valued function of k whose jump discontinuities lie 
along a system, S, of contours including (generally) the real and imaginary fc-axes. Second, one introduces 
a so-called g-function (an innovation discovered in |13|). a scalar function of k analytic for fc S C \ E, and by 
multiplying the matrix- valued unknown on the right by the diagonal matrix eS[k)a3/E ^^^^ arrives at a second 
equivalent Riemann-Hilbert problem on the same system S of contours but for which the jump conditions 
involve the boundary values of g. The function g is then chosen to cast the jump conditions for the matrix 
unknown in as simple a form as possible. By the phrase "as simple a form as possible," we really mean that 
g is chosen so that the steepest descent method of Deift and Zhou [TS] applies, implying a set of conditions 
that one expects to uniquely determine g. A matrix-valued Riemann-Hilbert problem to which the steepest 
descent method applies is one for which the ratio of the boundary values of the unknown is a matrix taking 
two alternate forms in complementary systems of subintervals of E. The number of subintervals and their 
common endpoints depend on the remaining parameters of the problem, including x and t. 

The common endpoints fc = kj[x, t) of these subintervals have special significance as branching points of 
the Lax eigenfunctions (solutions v of the linear equation (|176p ; see the appendix) corresponding to the local 
asymptotic solution of the MNLS equation in a small neighborhood of {x,t). If there are just four of these 
endpoints, then one expects the local solution to have plane-wave form (j){x,t) = A{x,t)e^^^^'*^^'^ + 0(e), 
where one may consider p = and u ~ dS/dx to be constants. To understand the interpretation of the 
branching points kj in terms of p and u, note that the relevant scattering problem (see (j47[) below) is written 
in terms of p and u as 

-4fc2 + 1 - au Aap^/^k 
Aap^/'^k 4fc2 - 1 + Q 

With p and u constant, this is a constant-coefficient linear problem with cigcnsolutions w(.t) = e^±^/^w±, 
where 

-4fc2 + 1 - au ^ap^f^k 
Aap^^'^k 4fc2 - 1 + au 
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dw 
2a£— 
ax 



(22) 



w± = 2aE±^/v± . (23) 



The eigenvalues are 



2aE± = ±iy^l6a^pk^ + (4/fc2 + au - I f , 



(24) 



and branching points of the eigenfunctions in the complex fc-plane will occur where £'+ = or equivalently 
where E± = 0. There are therefore at most four branching points satisfying 

4fc| = -2aV -au + l± 2oLyJ{o?p + au - l)p . (25) 

It follows that if (i) steepest descent analysis of the inverse-spectral Riemann-Hilbert problem in the limit 
e I leads to a solution of the MNLS equation of the form cj) = A{x, f)e''^(^'*)/^ +0{e) where p{x, t) = and 
u{x,t) — dS/dx, and (ii) the plane wave characterized by p and u is modulationally stable according to the 
hyperbolicity criterion ([2T|) . then it is necessary that the ^-function of steepest descent theory provides four 
endpoints kj {x, t) all lying on the imaginary axis. On the other hand, if p and u characterize a modulationally 
unstable plane wave, then the four endpoints kj {x, t) will need to form a quartet invariant under the reflections 
k ^ k* and k — > —k*. 

The unstable case of course requires that the contour system, S, contain components that are disjoint 
from the real and imaginary /c-axes, and determining the correct boundary values of the function g on such 
a contour further requires comparing with certain spectral functions computed on the real and imaginary 
fc-axcs (see the appendix) that are analytically continued to the complex /c-plane. This process of analytic 
continuation, that we expect to only be required in the unstable case, is the spectral analogue of analytic 
continuation of the functions p(a;,0) and u{x,Q) to the complex x-plane as would be necessary to solve the 
Cauchy initial-value problem for the elliptic modulation equations. 



2.4. Other forms of the modulation equations. 



2.4.1. Local conservation laws. With the dependent variable v := pu + ap^ used instead of u, the system 
(Uni) takes the form 

, ^ ( „ 2 3a 2 



(26) 



dv d f V p- , 

dt + d-x[i;+^'''-Y)-'- 

That is, (|19p may be written as a system of local conservation laws for densities p and v. Strictly speaking, 
the system ([T9)) is already in this form: 



dp d f 3a 2 



du d f 1 



apu — p + -u = ; 



dt • dx V 2 



(27) 



however since u = dS/ dx, the local conservation law for u is trivial. Consequently one should seek a second 
conserved local density that is not a derivative with respect to x of a local expression in A and 5*. It 
is well-known and easy to verify that the momentum Vq := pu is a nontrivial conserved local density for 
the modulation equations (jlip associated to the focusing NLS equation, and the nontrivial conserved local 
density v for (|19p (or equivalently, (|26p) is clearly an a-perturbation of the momentum. 

In the paper [16j the quantity v is called a noncanonical momentum. The authors of that paper give 
an analysis of the (equivalent, by a Galilean boost with speed c = a~^) derivative NLS equation in the 
semiclassical limit. In the spirit of the paper of Grenier |20], the authors use the system ([26|) (with u replaced 
by M -|- as is appropriate for the derivative NLS equation rather than the MNLS equation) to model the 
behavior of ipix, t) before singularities form in the solution of the limiting system. The development in the 
paper [16| appears to assume, however, that the system of conservation laws (|26p is hyperbolic for arbitrary 
initial data, and we have seen (see (|2ip which becomes simply p + au > when u is replaced by u + a~^) 
that this is not the case. 



2.4.2. Connection with focusing NLS in the semiclassical limit. Suppose that a^p + au — 1 < and that 
p > 0, defining an open domain D^{ol) C for each a > 0. Then consider the map : {p,u) G Z3_(a) i-^ 
(p, u) G M."^ defined by the formulae 

p = —p ■ {c? p + — 1) , it = u + 2q;p . (28) 

The mapping is one-to-one and maps I?_(a) onto the upper half-plane 

R-^{ {p, u) e such that p > 0} . (29) 

The inverse FZ^ : i?_ D^(a) is given by the formulae 

p = ^ (aii - 1 + v/(aw- 1)2 -h4a2p) , u ^ ^ (l - {au ~ 1)^ + ia^ pj , (30) 
where the positive square root is meant. The mapping F- is illustrated in Figure [1] 

P 





Figure 1. Top: the domain of F^ is shaded. Bottom: the range of F^ is shaded. Note 
that Q{u; c, a) := [(c - a"'^)'^ - {u - a"i)2]/4. 

It is a direct calculation to check that when the MNLS modulation equations are restricted to functions 
{p{x,t),u{x,t)) £ D-{a), they take the form 
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in the new variables {p,u). Comparing with for /c = 1 we see that on the modulationally unstable 
sector of the phase space the formal semiclassical limit for the MNLS equation is, when viewed in the correct 
variables (p, it), the samcQ as that for the focusing NLS equation. 

■^Note however, that the mapping F— docs not appear to estabhsh any equivalence between the MNLS equation and the 
focusing NLS equation at the level of finite nonzero e. We do not know whether the e-independent mapping _F_ can be prolonged 
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This identification allows us to immediately apply certain known results for the focusing NLS equation to 
the MNLS equation. For example, a particular solution to (j31[) was obtained in implicit algebraic form by 
Akhmanov, Sukhorukov, and Khokhlov [5] in 1966: 

p = {Al + t'^ p^) sech^ {x - ut) , u = -2</5tanh(x - ut) . (32) 

Here > is an amplitude parameter, and of course x and t are the independent variables in (|3ip . Setting 
t — gives the initial conditions 

p{x, 0) = Al sech^ {x), u{x,0) = 0, (33) 

(which corresponds to a "chirp-free, return-to-zero" pulse in the optical fiber context) and for small t the 
Implicit Function Theorem applies to (j32p and allows us to solve for p{x, t) and u{x, t). The earliest positive 
t for which the Implicit Function Theorem fails is t = {2Ao)~^, and the singularity occurs for a; = 0. 
For the implicit solution of Akhmanov, Sukhorukov, and Khokhlov, p{x, t) is an even function of x and 
u{x,t) is an odd hmction of x for each t e (0, (2Ao)~^). For any a > it is possible to compose the 
Akhmanov-Sukhorokov-Khokholov solution with the transformation FZ^ and thus obtain a solution of the 
MNLS modulation equations that remains restricted to the modulationally unstable sector of the phase 
space. The transformation FZ^ breaks the even/odd symmetry of the functions p and u. Indeed, the initial 
conditions (1551) both become even midcr F^^: 



p{x, 0) = - 2^ (^1 - \/l + 4a^^o sech2(x) j , u{x, 0) = ^ ^1 - y^l Aa^Al sech^(x) j . (34) 

These have the following asymptotics in a: 

p{x,0) = Alsech^ix) + 0{a^) , u{x,0) = -2aAlsech^{x) + 0{a^) , as a ^ (35) 

and 

p{x, 0) = Aoa^^ sech{x) + 0{a^^) , u{x,0) = -2AoSech{x) + 0{a^^) , as a ^ oo , (36) 

which shows that the MNLS pulse corresponding to the Akhmanov-Sukhorukov-Khokhlov solution has "fatter 
tails" when a is larger, but is of lower amplitude. The even symmetry of p and u at i = is, however, 
irrelevant under evolution. To show the dynamics of the Akhmanov-Sukhorukov-Khokhlov solution composed 
with FZ^ for different a we have numerically inverted the implicit relations (j32p and composed explicitly 
with FZ^, plotting the solution snapshots for six equally spaced times tk G [0, (2^o)~^] for Ao = 2 and 
a = 1/2 (Figure [2]) and for Aq = 2 and a = 2 (Figure [3]). In both cases we see leftward-leaning pulses 
evolving toward a singularity at t = {2Ao)^^ = 1/4 ^md x ~ 0. 

2.4.3. Connection with defocusing NLS in the semiclassical limit. Riemann invariants. Now suppose that 
a^p + au — 1 > and p > 0, defining an open domain D^{a) C for each a > 0, and consider the map 
F^ : {p,u) e D+{a) t-^ {p,u) G defined by the formulae 

p ~ p ■ [o? p + au — 1) , u — u + 2ap . (37) 

The range of the mapping F^ is a strict subset of the upper half-plane: 

R+{a) = {{p,u) e such that < p < (u - a^^f/A}. (38) 



into a gauge transformation relating the MNLS equation to the focusing NLS equation for general e. However, in [22] (see also 
I35p Hayashi has shown that if r) is a solution of the derivative NLS equation 

■-— H ^ -I- tae — (1,/.!^. 

dr 2 di 

(which is itself equivalent to the MNLS equation via a Galilean boost with velocity c = 1/a as discussed in the introduction) 
then 



ation via a Galilean boost with velocil 
/:=V^cxp(^H^ 1^ m^',T)\^di'^^ and c? := v^exp (^^^ l^(«',r)|2<) 



9? 



satisfy the coupled system 



^^^^-1-19 =0 and + — ^ttk + W J 



' dr ' 2 dr 2 

which reduces to the focusing NLS equation if g = —if and to the defocusing NLS equation if 3 = if ■ However neither of these 
reductions is consistent with the definition of / and g in terms of 1/1. 
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Figure 2. Snapshots of the Akhmanov-Sukhorukov-Khokhlov solution of the focusing NLS 
modulation equations interpreted via FZ^ as a solution of the MNLS modulation equations 
for a ~ 1/2. The initial conditions are shown with blue curves and the snapshot correspond- 
ing to the singularity formation time t ~ \/A is drawn with red curves. 




p' 


+ 


u 


p 


d 


P 


u 




1 


u 


dx 


u 



0, 



(39) 



Unlike the map is not one-to-one. Each point in the interior ol R^(a) has exactly two preimages 
in _D_|_(a), one with u < oT^ and one with u > a~^. On the other hand, each point on the boundary curve 
p = (u — a"^)^/4 has exactly one preimage in _D+(a), satisfying u = a^^ . See Figured! 

Again, by direct calculation, the MNLS modulation equations (|19p for p{x,t) and u{x,t) are transformed 
under F^ into the system 

d_ 

dt 

which upon comparison with (jlip for k = — 1 we recognize as the modulation equations for the defocusing 
NLS equation. Once again, we are led to expeclQ that in the semiclassical limit of e 4 0, the dynamics of the 
MNLS equation on the stable sector of its phase space should be equivalent (under the mapping F+ ) to the 
dynamics of the defocusing NLS equation. 

In principle, every hyperbolic quasilinear system with two dependent variables can be written in diagonal 
form through the introduction of new dependent variables called Riemann invariants. A suitable definition 
(see, for example, [26]) of Riemann invariants for the system ((39|) is 



(40) 



with the same caveat as indicated in the previous footnote 



10 



-2 -1 





X 



1 2 




-1 
-2 

u 

-3 
-4 



-2-10 1 2 

X 

Figure 3. Snapshots of the Akhmanov-Sukhorukov-Khokhlov solution of the focusing NLS 
modulation equations interpreted via FZ^ as a solution of the MNLS modulation equations 
for a — 2. The color scheme is the same as in Figure\^ 




where the positive square root is meant. Indeed, in terms of these variables, the system (p9| takes diagonal 
form: 



dR+ 
~dt 



The quantities 



+ 



R+ 
1 

~ 2 



R- 



3R- 



dR+ 

dx 
dR_ 

dx 



0, 



0. 



3i?± + R^ 



(41) 



(42) 



are the characteristic velocities^ that is, the eigenvalues of the coefficient matrix of (|39p written in terms 
of the Riemann invariants. The parabolae in the range i?+(a) of as shown in the bottom diagram of 
FigurelH are coordinate curves for the Riemann invariants; indeed, the relation p = [u — c)^/4 corresponds 
to the ray i?+ = c/2, R^ < c/2 if u < c and the ray i?_ = c/2, i?+ > c/2 if u > c. 

Diagonal form made possible by the introduction of Riemann invariants is useful for many purposes. For 
example, the way wc discovered that the transformation transforms the MNLS modulation equations 
((T^ on the modulationally stable sector into the defocusing NLS modulation equations ([55]) was by first 
writing (I19[) in diagonal form and subsequently comparing with the known diagonal form of (j39p . (Later we 
generalized this result to obtain the corresponding result that F_ links (|19[) on its modulationally unstable 
sector to pip although as both are elliptic systems real Riemann invariants cannot exist as an intermediate 
step.) Another advantage of diagonal form is that one can easily see that there are nontrivial solutions for 
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Figure 4. Top: the domain of is shaded. Bottom: the range of is shaded. 

which either i?+ or i?„ is identically constant. Such solutions are called simple waves and their dynamics 
are characterized by the scalar inviscid Burgers equation. Indeed, if R± ~ c/2, then it is easy to show that 
(Tq: evaluated for R± = c/2 and i?zp arbitrary satisfies 

d(T± da± 

whose general solution with arbitrary initial condition (t±(.t,0) = f{x) can be written in implicit form as 

'j± = fi^ - <y±t) ■ (44) 

As is well-known, if / has critical points, then this solution only exists in the classical sense for a finite 
time due to the formation of shocks (derivative singularities). The number of shocks that can form is 
bounded above by the number of critical points of / (whether or not more than one shock forms is a 
sensitive question of regularization that we do not address here; in particular we arc ultimately interested in 
dispersive regularization where an early-occurring shock can "spread" into microscopic oscillations and wipe 
out other shocks that might occur in nonconservative regularizations). A final advantage of diagonal form 
that is useful in considering solutions more general than simple waves is that as long as a solution remains 
classical it is bounded in terms of its initial data by sharp inequalities of the form 

inf i?±(y,0) < R±{x,t) < supi?±(y,0). (45) 

In other words, the smallest coordinate box containing the solution at a given t > is in fact independent 
oft. 

Perhaps with the help of Riemann invariants, the mapping F^ and its double- valued inverse F^^ can 
be used to obtain solutions of (fT9|) from known solutions of the hyperbolic system (|39|) . a problem that is 



better understood than the mixed- type system p^ . having received much attention in the hterature (see, 
for example, [4]). As the mapping is not one-to-one, some care must be taken in using given initial data 
for in its modulationally stable sector to obtain corresponding initial data for ([55)1 and then porting the 
solution thereof back to the domain D+(a) via F^^. 

Consider, for example, smooth pulse-like initial data for for which p{x, 0) ^ as x ^ ±oo and for 
which the modulational stability condition (pij) holds for each x £ E. If additionally u{x,Q) > a^^ holds 
for all x e R, then the image of this initial data in the (p, u) plane is a curve parametrized by a; G K. that 
is attached to the p = axis at one or more points with u > and that avoids the boundary curve 
p ~ {u ~ q;~^)^/4 of the range i?+(a). As parabolic curves of the form p = {u — c)^/4 for c > are 
coordinate curves of the Riemann invariants for p9p , the evolution in time prior to shock formation will not 
alter this basic picture due to the inequalities (|45p . Thus, although when applying the inverse mapping F^^ 
pointwisc in x one has to make a choice between two distinct preimages, there are only two distinct pairs 
of continuous functions {p{x,t),u{x,t)) for each t corresponding to the evolved initial data. Furthermore, 
exactly one of these will satisfy the boundary condition that p{x,t) — > as a; — *■ ±cxd; the remaining solution 
of (|19p will satisfy a^p{x,t) -\- au{x,t) — 1 — > as x — > ±oo instead, and taking into account the initial 
condition this extraneous solutior[j is not continuous in time t. On the other hand, if for some finite x e M 
the initial data crosses the threshold into the subregion of the domain D+{a) for which u < a~^, then even 
the image of the initial data as a curve in the (p, u)-plane may have some unexpected properties. Indeed, 
a direct calculation shows that all curves in the (p, u)-plane passing through the vertical axis u = 
with slope dp/du ^ —{2a)~^ will have images in the (p, u)-plane that are tangent to the boundary curve 
p ~ {it — a~^)'^/4: while this is not the case for curves passing through with dp/du — — (2a)~^ (as holds for 
all of the coordinate parabolae). Therefore, if dp/du = — (2q!)^^ for u = a^^ then the transformed initial 
data functions p(x, 0) and m(x, 0) will necessarily have a common critical point for some x € M. Moreover, 
for initial data crossing the u — a^^ threshold, there exist multiple preimage function pairs {p{x,t),u{x,t)) 
that are continuous in x with the same boundary conditions as x ^ ±cx3. Care must be taken to select the 
correct preimage to have continuity for t > 0. 

These facts suggest that while the modulationally stable sector of the MNLS modulation equations (fT9|) 
can be mapped to the defocusing NLS modulation equations ([39]) which we may consider to be a problem for 
which some intuition is available, new and important dynamical features of the MNLS modulation equations 
can be introduced simply through the change of variables F+. For example, let c > and j4o > be 
fixed, and consider initial data for (flQl) of the form 

p{x,0) = Alsech^{x) , u{x, 0) = c - 2Ao\/ac ~ 1 • sech(x) . (46) 

The corresponding curve parametrized by x G K lies along one of the coordinate parabolae in the (p, u)- 
plane with u < c and p — > as x ^ ±oo. If < < (c — a~^)/{Aa), then the initial data satisfies 
u > l/a for all x and thus a single-peaked pulse is mapped to a single-peaked pulse as shown in Figure [S] If 
(c— a~^)/{4a) < < {c — a~^)/a, then the initial data crosses the threshold of u = but does not yet 
reach the point of tangency between the coordinate parabola and the stability boundary a^p + au —1 = 
(which is also mapped to p = along with p = 0). In this situation, the image in the (p, ii)-planc of the part 
of the initial data with u < is the same as that of part of the initial data with u > a~^, and this results 
in a double-peaked image pulse as shown in Figure [H] Finally, if exceeds (c — a^^)/a, then the part of 
the coordinate parabola in the (p, {t)-plane that is retraced is the entire left-hand branch for u < c, and part 
of the right-hand branch (along which p can become arbitrarily large) is also traced out in the center of 
the pulse. Therefore, in this situation the single-pulse initial data for the MNLS modulation equations 
is mapped by F^ to triple-pulse initial data for the defocusing NLS modulation equations ((39)) . as shown 
in Figure [71 Since the coordinate curves correspond to constant values of one or the other of the Riemann 
invariants, the dynamics of all of these initial conditions under the defocusing NLS modulation equations 
P9p will be given by simple waves. It follows that the evolution is governed by the scalar inviscid Burgers 
equation, and it may thus be shown that the number of shocks that can form corresponds to the number 
of critical points of the simple- wave initial data. Thus, simply by varying the amplitude of a single-peaked 
pulse- like initial condition for the MNLS modulation equations (|19|) , one can produce dynamics that lead to 
either one, two, or three (dispersively regularized) shocks. 



'Of course this extraneous solution may be meaningful if different boundary conditions are imposed at x = itoo. 
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Figure 5. The initial data and its image under F+ for a = 1 and c = 6 with = 5/8. 
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Figure 6. The initial data and its image under for a = 1 and c = 6 with = 30/8. 



A relation similar to the equivalence between the MNLS modulation equations (|19p and those for the 
focusing or defocusing NLS equations has been observed by Kuvshinov and Lakhin |30| . Working with 
the equivalent derivative NLS equation, Kuvshinov and Lakhin derive the modulation equations (Whitham 
equations) corresponding not to complex exponential plane waves (genus zero) but rather to periodic waves 
given by elliptic functions (genus one). By comparing the Riemann invariants for these Whitham equations 
with the corresponding well-known formulae for the NLS genus one Whitham equations they deduce complete 
equivalence between these two systems of modulation equations. No doubt this phenomenon extends to 
genera greater than one, and a deeper explanation of these coincidences surely exists at the level of the 
algebro-geomctric description of arbitrary genus multiphase wave solutions of the NLS and derivative NLS 
equations. We will treat this problem in this framework in a future publication. 
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Figure 7. The initial data and its image under Fj^ for a = 1 and c = 6 with 



3. Bounds on the Discrete Spectrum 

In this section we establish some basic conditions on the discrete spectrum associated to the spectral 
problem ()176p for the MNLS equation. As seen in the appendix, for each k with Im{fe^} ^ there are two 
one-dimensional subspaces of solutions that decay to zero as x ^ ±(X) respectively. We refer to eigenvalues 
as values of k for which these two subspaces coincide and eigenfunctions as the corresponding solutions v 
of (|176|) associated to each eigenvalue that decay to zero in both direction^. Eigenvalues are exactly the 
zeros of of the function 5'ii(fc) in the region Im{fc^} < and thet zeros of the function S22{k) in the region 
Im{fc^} > 0, where S'ii(fc) and 5*22 (fc) are defined by (jl84p . The integrabililty of the MNLS equation leads to 
the fact that eigenvalues do not vary in time (see (|25ip ). Consequently, we consider the spectral problem at 
t = Q and set as the initial condition for the Cauchy problem for the MNLS equation ([T]). Additionally, we 
note that due to the symmetries (|230p and (|231[) , the eigenvalues have a four-fold symmetry, where if any of 
fc, fc*, — fc* or —k is an eigenvalue the others are as well. This symmetry allows us to restrict our attention to 
(say) the first quadrant of the complex fc-plane. Any result established for k in the first quadrant therefore 
immediately gives analogous results for k in any other quadrant. 

We will be considering the linear differential equation (|176|) with a potential of the form (l>{x) = A(a;)e*'^^^^/^ 
where A{-) and S{-) are real- valued functions such that A> Q and A and S" decay rapidly to zero for large 
At first blush, one difficulty that presents itself is that there are rapid oscillations in the coefhcients 
of p76|) when e is small. This difficulty is, however, circumvented by the following device. If v is a vector 
solution of (fT76)) . then ^^ := e-^^'^^^^/^^^^'^v satisfies: 

2ae^ = iMw , (47) 
ax 



M 



(48) 



±oo 



where 

" -4fc2 + l-a5"(a;) AakA{x) 

4aA:A(x) ^k^-l + aS'{x) 

Since we have only scaled v by exponential factors with purely imaginary exponents, w decays as x 
if and only if v docs, that is, the discrete spectrum of (|47)) is the same as that of (|176p . 

A formal approach to (j47|) when e is small is apply the WKB method; that is, to assume a solution w 
of the form w = e*'^/'-^"^-' (wq + ewi + •••). Substituting this expansion for w into ([17|) and equating terms 
of the same powers in e gives a hierarchy of equations relating a and {w/;}. The leading order (eikonal) 



'^This is inexact usage, as tlie spectral parameter k appears nonlinearly in the system I I176I I. The concept is, however, 
completely analogous to the usual situation. 
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equation in e is: 

Mwo = , (49) 

and, consequently, da/dx must be an eigenvalue of M. This leads us to consider the eigenvalues of M as 
smooth functions of x. 

A simple calculation shows that the eigenvalues of M are ±w where 

uj[x- k) := [l&a'eAixf + {Ak^ - 1 + aS'{x)f] . (50) 

In WKB theory, the so-called turning points are those values of x for which the eigenvalues coincide: da/dx = 
±w = 0. Near such points, WKB theory breaks down as infinite derivatives appear and approximate solutions 
become multiply-defined. Setting a;(x; k) — 0, we see that for a given fc S C the turning points x G M are 
the solutions of the equation 

'iiakA{x) = ± (4A:2 _ l + aS' {x)) . (51) 
On the other hand, as a; S M varies, the equation (|5ip defines a parametrized curve T in the four quadrants 
of the complex fc-plane that we will call the turning point curve associated with the functions A{-) and S'{-). 
By definition, fc € T if and only if Imjfc^} ^ and there exists at least one real turning point x. As ((5T|) is 
quadratic in k, we may solve for k and therefore explicitly present T as a parametrized curve: 

k{x) = ^ {siaA{x) + S2\/a^A{xY +aS'{x) - l) , (52) 

where Si and S2 are independent signs (sj = ±1) yielding up to four branches of T . 

At this juncture we may make the observation that if a; is a point at which the hyperbolicity (modulational 
stability) condition (PT|) holds (recall that p := A? and u := 5'), or in the degenerate case if a^p + au— 1 = 0, 
then the four fc-values produced by the formula (|52p are purely imaginary and by definition do not lie in T. 
On the other hand, if x is such that a'^p -f au — 1 < 0, then ([5^ gives four distinct points k 



Im{A:} = si—A{x) . 



Re{fc} = S2-^Jl-aS'{x)-a^^A{xf 



(53) 



The equations ([55]) when restricted to the values of x for which the strict inequality a^p -I- au — 1 < holds 
(x-values of modulational instability) thus give a parametric representation of the turning point curve T. In 
particular, if the nonstrict modulational stability condition a^Aix)^ + aS'{x) — 1 > holds for all x S R, 
then by definition the turning point curve is empty: T = 0. 

Throughout this section, we will restrict our attention entirely to those values of fc S C with Imjfc^} ^ 
that do not lie on the turning point curve. For such fc, the function u}{x\ k) is well-defined by ([50)1 if we insist 
on continuity for all x £ M and impose the boundary condition that 

lim uj{x;k) =Ak^ -l + aS', where S' := lim S'{x). (54) 

Suppose that A{-) > is of class L^(IR) and is uniformly Lipschitz. In particular, this implies that A{x) decays 
to zerc3 as X ^ ±oo. Suppose also that 5"(-) is uniformly Lipschitz and that 5"'(-) is of class L^(R) (assuring 
the existence of limiting values of S'{x) as x ±oo). The function uj{x; k) defined as described above then 



''It is enough to consider the behavior as a; — > -|-oo. Arguing the contrapositive, suppose that A{x) > is integrable on 
R and Lipschitz, and that hmsup^^_|_^ j4(x) = M > (wc admit the possibiUty that M = +oo). We will show that the 
Lipschitz constant of A is bounded below by an arbitrarily large number, and hence is infinite. Begin by choosing 5 > to be 
arbitrarily small. Let {xn\ be a sequence of real numbers tending to -|-oo such that limn^oo A{xn) = M > 0. By passing to a 
subsequence if necessary, we may assume that Xn+i > Xn + & for all n. We first claim that In '■= ^'^fx„<x<x„+S ^{^) tends to 
zero as n ^ oo. Indeed, were this not the case there would exist some e > and some rio such that In> £ holds for all n > no- 
On the one hand we would then have 

A{x) dx > eS , n > rio , 



f 



and on the other 



/ A(x)dx< / A{x)dx<oo, 



which yields a contradiction. Since A is continuous, there exists for each n a number y„ S [x„,x„ + 5] such that = /„. 

As /„ tends to zero as n — > oo, for all n sufficiently large wc will then have both A(x„) > 2A//3 (or > 2 if M = -|-oo) 
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has the following properties whenever Im{fc^} 7^ and k ^ T: Lo{x;k) and its a;-derivative uj'{x;k) are 
uniformly bounded, and io{x] k) is uniformly bounded away from zero. Indeed, given the assumed properties 
of A{-) and S'{-) it suffices to show that uj{x; k) is bounded away from zero. But this is so since for ^ T 
we have uj{x; k) ^ for all a; G M, while the boundary condition and the limit 

lim uj{x;k)=±{Ak^ -l + aS'_), where S'_ := lim S'{x), (55) 



X — > — oc> 



when combined with the condition that Im{fc^} 7^ shows that u!(x;k) has nonzero limiting values as 
X ±00. 

Given uj{x; k) defined as above, we set 

2akA{x) d ( Mpc) \ 

cuix;k) 'd^:''^[io{x;k)+4k^-l + aS'{x))- 

Supposing the same conditions on A{-) and S'{-) as above, the function q{x;k) is uniformly bounded for 
X G M whenever Im{fc^} 7^ and k ^ T and a > 0. Indeed, writing q{x; k) in the equivalent form 

( i.\ J^fA'(\ uj'{x-k)+aS"{x) \ 

q{x;k)= [A{x)-'A{x) .,2 — 1 . Qu \ l ' 

uj(x;k) \ uj(x;k)+4:k'' — l + ab'(x)J 

our attention falls on the denominator in the second term: D := ll){x; fc) + 4fc^ — l + aS'{x). This denominator 
cannot vanish for any x G M because this would imply that uj{x; k)^ = (4fc^ — 1 + aS'{x))'^ which by (jSO]) is 
in contradiction with the assumptions that Im{fc^} 7^ 0, A{x) > 0, and a > 0. According to the boundary 
condition ([5^ . D also has a nonzero limit as x +00 for Im{fc^} 7^ 0. A similar situation holds as a; ^ —00 
as long as the limiting value of ti;(x; k) is achieved with the "+" sign in ([SS]) . It follows that q{x] k) is bounded 
for large |x|, at least if the limit ([55]) holds with the "+" sign. If (|55p holds with the "— " sign, then we must 
expand D as x ^ —00: 



whence it follows for Im{fc } 7^ and a > that 

so that in this limit q[x; k) is again bounded for large negative x. 

Finally, in terms of uj{x] k) and q(x] k) we define, for Im{fc^} 7^ and k ^ T, the quantity 



sup 



1 



2 sup 



Rc{q{x-k)} 



Im{ti;(a;; fc)} 



(60) 



dx \Im{a;(x; fc)}^ 
It is possible to have ~ +00. 

3.1. The fundamental condition on eigenvalues. Our main result is the following. 

Theorem 1. Let A : M. ^ M_|- he a uniformly Lipschitz function of class L^{M.) and let S' : M. ^ M. be 
uniformly Lipschitz with S"{-) of class L^(U.). Let k be a fixed complex number with /m{fc^} 7^ and k ^ T. 
If (j){x) := ^(a;)e''^'^^-'/^ is the potential in the linear differential equation (|176p . then the following statements 
hold: 



(a) // k is an eigenvalue, then 

(b) If Lk < {ae)~^ then k is not an eigenvalue 



|Im{/fc}| < ^sup^(x). (61) 

2 xeR 



and A{y„) < A//3 (or A(j/„) < 1 if M = +00) while \x„ — y„\ < 5. Therefore, for all such n, 

\A{xn) - A(y„)\ ^ Cm_ 

\Xn-yn\ ~ <5 ' 

where Cm is a finite positive number depending only on M . As & was arbitrarily small, this shows that the Lipschitz constant 
of A is arbitrarily large, and hence infinite. 
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Our proof of this theorem was inspired by some notes [T2] we obtained from Percy Deift regarding his proof 
with Stephanos Venakides and Xin Zhou of the corresponding "shadow bound" estimates for eigenvalues of 
the nonself adjoint Zakharov-Shabat problem. (We will explain more about their result, and show how it can 
be deduced from ours as a special case, at the end of this section.) 

Proof. To prove part (a), suppose that k is an eigenvalue and w is the corresponding eigenfunction. Using 
(|Tf)) . we calculate 

2ae (^wVa^^ = -i {Ak^ - 1 + aS'{x)) |wp - 8akA(x) Im{?i;>2} , (62) 

where |wp := + 1^2^, and we adopt the notation that will be used throughout the remainder of the 
paper that wi and W2 are the first and second components, respectively, of the vector w. As shown in the 
appendix (see (|198|) and (|217p ). the condition that A{-) £ L^{R) guarantees that the eigenfunction w 7^ 
decays exponentially to the zero vector as a; ^ ±00. Integration by parts then establishes the following 
identity: 

/ w^a^—j—dx^ / [w[wl — W2W2] dx = — i / [w[wl — w 2102]* dx ] , (63) 
J —00 "'•^ J —00 \j —00 / 

so that all three expressions are purely imaginary. Therefore, integration of (j62p and taking the real part 
yields 

/+00 /•+00 
|w|^dx- aRe{fc} / A(x)\ni{wlw2} dx . (64) 
-00 'J —00 

Since Re{fc} 7^ 0, jH gives: 

/+00 
A{x) Im{wlw2} dx 

inii/^| = ^25L__ . . (65) 

|w| dx 

> 

Applying the inequality 2|Im{w*W2}| < 2|i/;i||w2| < + |^^2|^ ~ |w|^ we obtain: 

|Im{A:}| < ^sup^(a;), (66) 

which proves part (a) of Theorem [1] 

To prove part (b), suppose that k is an eigenvalue of (|176p . or equivalcntly, of (|T7|) . We will now suppose 
that ifc < (ae)~^ and derive a contradiction. 

Let T(a;; fc) be a 2 x 2 matrix with eigenvectors of M corresponding to (distinct) eigenvalues uj{x; k) and 
~Ld{x\ k) as its first and second columns respectively. Note that T(a;; k) is invertible for all .t e M because k 
does not lie on the turning point curve T. Make the gauge transformation p = T(x; fc)~^w, where by means 
of ((17|) . p satisfies 

2ae-^ — — 2aeT(x; fc)^"^— — (x; fc)p + itoix; k)a3p . (67) 
dx ax 

A general eigenvector matrix is only determined up to multiplication on the right by a diagonal matrix, and 
we fix the gauge so that T(x; k)~^ ■ T'{x; k) is an off-diagonal skew-symmetric complex matrix depending 
on A{x), S'{x) and k. Thus, with the particular eigenvector matrix 

4:akA{x) H{x;k)' 
H{x;k) ~2uj{x;k) 

H{x; k) AakA{x) 
2Lo{x;k) H{x;k) 

where H{x; k) is a continuous function of x satisfying 



T{x-k) 



(68) 



H[x; kf = 2w(x; k) {uj{x\ k) + Ak^ - I + aS'{x)) , (69) 
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we see that T(x; k) ^T'(.t; k) is off-diagonal as desired: 

dT 

T{x;k)-^ — {x]k)^iq{x;k)a2, (70) 
ax 

where q{x; k) is defined by ([55)) . Moreover, det(T(x; k)) = 1, and, by the same arguments as were used to 
deduce the boundedness of q{x;k), the elements of T(a;; fc) are uniformly bounded in x (and hence so are 
those of T(x; k)~^). Therefore, (p7)) becomes 

2ae^ = iuj(x; k)(73p — 2iaeq(x; k)<72P ■ (71) 
ax 

Next we differentiate pVap with the help of (j7ip and the identities cr| = I, a2(J3 = i(Ji, and (T3(T2 = — ici 
to find: 

2ae-^ (p^o-3p) = -2 Im{cj(a;; fc)} Ipl^ - 4ae Rc{q(x; fc)}pVip . (72) 

ax 

Since T(x;fc)~^ is uniformly bounded in x, p decays to zero for large |a;|. Dividing through in (|72p by 
Im{a;(a;; k)} and integrating gives 



" , ,2 , as d ^ /■+°° 2aeRc{q{x-k)} ^ 

, 'P' '^^-1-^ lmMx,k)}Tx^^'^^^^''-L M.(.;fc)} P'^^P^^ 
/■+°° _d_ / \ _ /■+°° 2a£Re{g(x;fc)} 

7-00 da; VIm{a;(x;fc)}/ P ^'^P ^ Im{a;(a;;fc)} ^ 

where the second equality is due to integration by parts and the decay of p. Since 



(73) 



Ip^ctspI 



bi|'-b2|' <|p|' and \p^<Jip\^\p;pi+pIp2\<2\p,\\p2\<\p\\ (74) 



we therefore obtain 

\p\'^ dx < aeLk / |p|^fia;, (75) 

where Lk is the constant, independent of e, defined in (pOjl . Cancelling the nonzero L^(R)-norm of p, we 
arrive at a contradiction. Therefore k cannot be an eigenvalue. This establishes part (b) and completes the 
proof of Theorem [TJ □ 

3.2. The shadow region. For A{-) and S{-) satisfying the conditions of Theorem [U wc have seen that 
uj'(x]k) and q{x;k) are uniformly bounded functions of x G M whenever fc is a complex number with 
Im{fc^} 7^ and k ^ T. It follows easily that under these conditions, for each such complex number k there 
exists a constant dk > 0, independent of e, such that 

Lk < dk sup [lm{uj{x; fc)}]~^ . (76) 

It is then of interest to ask for which values of fc can one conclude that Im{a;(a;; fc)} is bounded away from 
zero for all x e R, because Theorem [1] implies that such a value of fc cannot be an eigenvalue for any value 
of £ > that is sufficiently small. Note that Im{fc^} ^ implies that |Im{a;(x; fc)}| 4|Im{fc^}| ^ as 
X ±oo. Consequently, it is enough to show that Im{ti;(a;;fc)} 7^ for all a; e K. Wc will now derive a 
simple geometric condition in the complex fc-plane that implies that Im{u;(a;; fc)} is nonzero for all real x. 
For C € M, let Sc denote the hyperbola in the complex fc-plane defined by the relation 

Ec : {fc : Re{fc2} +C = Rc{kf - lm{kf + C = 0} . (77) 

Lemma 1. Given a fixed complex value ko in the open first quadrant of the complex plane, define the value 
Co — (Rejfco}^ — Imjfco}^). Let Iq he the set of x values that correspond to intersection points of the 
turning point curve T and Sp^. //2Im{fco} — aA{x) > for all x G do, then Im{[j(.T; fco)} 7^ for all x G M. 

Proof. For any complex number z, the condition Im{z} 7^ can be rewritten as a condition on z^: 

Im{z2} ^ 

Im{z} 7^ ^ <( or ) . (78) 

Im{z2} = and Rejz^} < 
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From ([50)1 wc find that lm{Lo{x; ko)^} is given by 

Im{a'(a;; hf} = 16 Im{fco}Re{fco} {2a^A{xf + aS'{x) - 1 - 4Co) • (79) 

Since ko is in the open first quadrant, Rc{fco}Ini{fco} > 0. Let Jo denote the set of a; G R sucli that 
Ini{w(x;fco)2} = 0. 

Case 1: Let a; e M \ Jq. Then by ([75)) we clearly have Im{LLi{x; ko)} 7^ 0. 

Case 2: Let x E Jo Consider the hyperbola Ec'o defined in ([77)1 . Each intersection point k of S^*^ with 
the turning point curve T corresponds to a real value of x for which a^A{xY + aS'{x) — 1 < according 
to the natural parametrization of the latter curve, given by (|53p . By eliminating Re{fc} and Im{fc} between 
the two equations of ([53)) and the definition ([77)) of the hyperbola Scoi it is easy to verify that each real 
X corresponding to such an intersection point k is an element of Jo- There may however be some x £ Jo 
that do not correspond to intersection points of with the turning point curve T. These are values of 
x for which k = k{x) defined by ([52)) lies on the hyperbola Ec„ but is a purely imaginary number because 
a'^A{x)'^ + aS'{x) — 1 > 0, and hence is not in T by definition. Since Im{fc^} ^ but lm{oj{x;ko)^} = 0, 
([79)) shows that such points x satisfy the completely equivalent inequality a^A{x)'^ — ACo < 0. Therefore, 
we define /q to be the subset of Jo consisting of a; S R such that 2a^A(a:)^ + aS'{x) — 1 — ACo = and 
a^A{x)'^ — ACo > 0. In other words lo consists of all x G Jo such that x also corresponds to a point of 
intersection of the turning point curve T and the hyperbola Sco- It is also clear, both directly from the 
inequality a'^A{x)'^ — ACo > and also from the fact that hyperbolae Sc with C < contain no purely 
imaginary points fc, that Co < implies lo = Jo- 
Let X be a fixed element of Jo and consider Re{a;(a:; k)^} restricted to fc-values lying on the hyperbola 
SCo- This restriction makes Re{a;(a;; fc)^} a function of Im{fc}^: 

Re{uj{x; kf} = -64Im{fc}^ + 6AColm{kf - Aa^A{xf{ACo - a'^A{xf) 

^•^^co (80) 

= -64 (Im{fc}2 - ri) {lm{kY - rj) , 

where 

n := ^a^A{xf , and r2 := ^{ACo - a^A{xf) . (81) 

As a function of Im{fc}^, the right-hand side of ([80)) is a concave-down quadratic. Case 2 now breaks into 
two subcases depending on whether or not x € Jo satisfies x £ lo- 

Subcase 2(a): Let x E Jo\Io- Assume that Re{a;(x; fco)^} > 0. Since obviously ko lies on the hyperbola 
Sc^, it follows from ([50)) that Imj/co}^ must lie in the closed interval whose endpoints are the two roots ri 
and r2, and in particular 

Im{fco}^ < max(ri,r2). (82) 

Eliminating Im{fco}^ using the definition of Co, we therefore find 

Re{fco}2 < niax(ri - C^, ^2 - Co) - (83) 

Now, r2 — Co < by definition. The assumption that x ^ lo however implies further that ri — Co < 0. There- 
fore, Re{fco} = which contradicts the assumption that ko lies in the open first quadrant. Consequently, for 
all X G Jo\Io, Re{uj{x] ko)^} < and therefore lm{uj{x; ko)} 7^ by ([75)) . 

Subcase 2(b): Let x G lo- It is always the case that ri > and for x G lo wc also have r2 < 0, by definition. 
Since (|80p is quadratic in Im{fc}^ and concave down, and since ko € ^c^, the inequality Imjfco}^ > ri implies 
that Ke{uj{x; ko)^} < for this value of x- Since x G lo C Jo, it follows from ([78)) that (for ko in the open 
first quadrant) 2Im{fco} — aA{x) > implies that lm{uj{x; ko)} 7^ 0. 

We have therefore shown that if a; e M \ /o, then Im{(jj(x; fco)} 7^ 0, while (according to subcase 2(b) 
above) if x G lo, we have Im{ti;(a;; ko)} 7^ as long as 2Im{fco} — aA{x) > 0, which completes the proof of 
the Lemma. □ 
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We might remark that Lemma [T] combines with the estimate ([75| and part (b) of Theorem [T] to give 
a weaker (e-dependent) version of part (a) of Theorem [T] because if 2|Im{fco}| > asnp^^^A{x) then the 
incquahty 2|Im{fco}| > aA(x) holds in particular for all x G Iq- 

Now define 

Cmin := i min (2a^A(x? + aS'(x) - l) C^ax :== \ max (2a^A(x? + aS'(x) - l) . (84) 

Corollary 1. Suppose A{-) and S'{-) satisfy the conditions of Theorem]^ and let ko be a fixed complex 
number such that Imlfc^} 7^ 0. If Co < Cmin or Co > Cmax then ko is not an eigenvalue for any e > 
sufficiently small. 

Proof. By means of ([79|) it is clear that if x G Jo, then 2a'^A{x)'^ + aS'{x) — 1 = 4Co. Therefore, either 
of the inequalities Co < Cmin or Co > Cmax implies that Jo is empty and thus (by Case 1 in the proof of 
Lemma[T]) we have lm{Lu{x; ko)} ^ for all .t G M. □ 

Note that Sc„i„ and Sc„ax are the osculating hyperbolae of the image of k{x) defined by (|52|) . which in 
turn is the union of the turning point curve T and (possibly) an interval X of the imaginary axis. This gives 
Corollary [1] geometrical meaning since the condition Cmin < C < Cmax is equivalent to the statement that 

i]cn(rux) 9^0. 

A refinement of Corollary [1] is an immediate consequence of Lemma [T] 

Corollary 2. Suppose A{-) and S'{-) satisfy the conditions of Theorem]^ and let ko be a fixed complex 
number such that l-ai{k^} ^ 0. //Cmin < Co < Cmax but Tic„ C\T = then ko is not an eigenvalue for any 
e > sufficiently small. 

Proof. The set lo is empty, so the proof follows from the estimate ((76|) and Lemma [TJ □ 

We may now formulate a more precise geometric condition to exclude the existence of eigenvalues. For 
each C G M with Sc n T 7^ we define a number Uc as follows: 

Uc sup Im{fc}2 . (85) 

Then set 

Sc := {k G Sc such that Im{fc}2 < Uc} , (86) 

and 

S:= U Sc. (87) 

C6[Cn,i„,Cn,ax] 

Corollary 3. Assume that A{-) and S'{-) satisfy the conditions of Theorem]^ Suppose that ko G C and 
Imlfcg} ^ 0. If also ko ^ S, then ko is not an eigenvalue for any e > sufficiently small. 

Proof. First note that, by Corollaries [1] and [21 it is sufficient to consider those ko ^ S for which Cmin < Co < 
Cmax and also Ec^ n T 7^ 0. Now, clearly ko lies on the hyperbola Eco- However, since Sc^ C S, we also 
have that ko ^ Sc„ , so 

Im{fco}2 > Uc„ ■ (88) 
By the definition of the set lo and the definition ([53]) of the turning point curve T, we see from (f85|) that 

Uc„ - sup . (89) 

Lemma [T] (and its obvious consequences in the remaining quadrants of the complex fc-planc) then implies 
that ko cannot be an eigenvalue as long as e > is sufficiently small. □ 



At this point we may observe an important consequence of the hyperbolicity condition (|2ip . 

Corollary 4. Suppose that A(-) and S{-) satisfy the conditions of Theorem[^ and also that o? A(x)^ + 
aS' (x) — 1 > holds for all x G M. The latter condition makes the modulation equations (|19[) (nonstrictly) 
hyperbolic for all x G M. at t — 0. If ko is any fixed complex number with Imlfc^} 7^ 0, then ko is not an 
eigenvalue for any £ > sufficiently small. 

21 



Proof. In this situation, the turning point curve T is empty, and hence so is the set S, so this statement is 
a direct consequence of CoroUary[3l □ 



The set S is the "hyperbohc shadow" of the turning point curve T. One can imagine hght being projected 
along the hyperbolae Sc from infinity, and the set S is the shadow cast by the turning point curve in this 
light. Let N be any open neighborhood containing the axes Im{A;^} = and the shadow S. Then from 
Corollary [3] and part (a) of Theorem [1] the set F := N D {k such that 2|Im{fc}| < asnp^^g^A{x)} contains 
all of the eigenvalues if e > is small enough. In particular, the set N may be taken to be the union of 
S with a thin "collar" surrounding the boundary of S and the axes. If one wants to contain the discrete 
spectrum for all < e < Emax, then the collar must be made thicker for larger Emax- The region F is perhaps 
best understood visually. Figure [8] contains two heuristic plots of F in relation to the turning point curve 
and the family of hyperbolae Ec- 




Re{k} 



Re{k} 



Figure 8. "Light'' is projected along the hyperbolae and casts the shadow on the turning 
point curve. The shadow region is shaded green. The collar is shaded grey and the red 
horizontal lines correspond to the threshold in part (a) of Theorem\^ 

The need for a "collar" arises from our methods of proof: wc found an upper bound for (see (|76p ) 
and argued that to ensure that the eigenvalue exclusion inequality < (ae)^^ holds for e small enough we 
simply checked for finitcness of the upper bound. The bound (j76p involved a constant dk that depended on 
k in an unspecified way. With better accounting of the fc-depcndence of upper bounds for Lk it is possible 
to obtain a sharper condition on eigenvalue exclusion that explicitly depends on e; in other words, one may 
estimate the width of the "collar" . As a simple example of such a calculation, we present the following result. 

Corollary 5. Assume the same conditions on A{-) and S'{-) as in Theorem]^ There exist constants > 
and Ca > independent of e, such that whenever \k\ > Ka and 



|Im{fc}| > 



|Re{fc}| 



(90) 



then k is not an eigenvalue. 



This Corollary says that the width of the collar about the real /c-axis scales like e^^^ when is large. 
Moreover, it shows that the collar becomes thinner when |Re{fc}| is larger. 

Proof. Another upper bound for Lk is clearly 



Lk < sup \uj'{x; k)\ ■ sup — -^-r-r^ 
xm 2;eK |Iin{w(a;; fc)}|^ 



2 sup \q{x;k)\- sup , . , — tttt 
xm X6B |Im{t^(a;; fc)}| 



(91) 



When k is large in magnitude, we have the following asymptotics holding uniformly with respect to a; G R: 
uj{x]k) = 4fc2 + 0(l), 

oj'{x; k) = \ {32a^ k^A{x) A' (x) + 2(4fc2 - 1 + aS' {x))aS" {x)) = 0(1) , (92) 
q{x;k) = 0{\k\-'). 

From the above it obviously follows that Im{w(a;; k)} = 4Im{fc^} + 0(1) uniformly in x for large k. When 
Im{fc^} is small, this estimate is not satisfactory, so to improve it we note that y/1 + x + iy = f + ig with 



/ = Wi(v/(TT^)^T^+(i + .T)) 



9 = sg: 



:n(2/)y ^ {,/{r+W+7 - (1 + x)) 



(93) 



where x and y are real and sufficiently small, and all square roots are positive. It follows that as z = x + iy 
tends to zero, 



VlT^ = I + Oi\z\) + iO{y) (94) 
where all order symbols denote real-valued functions. Applying this result to 



io{x; k) = 4fc^^l + (2«2^(x)2 + aS'ix) " 1)^ + (1 - ' (^5) 

and using the facts that 



we see that 

lm{uj{x; k)} = 4lm{k^} ■ (l + 0{\k\-^)) (97) 
holds uniformly for a; G M. We therefore have shown that 

Im{k'^}y \\k\ \Im{k^}\J \Im{k^}y \ \k\ Im{fc2}2 



Since part (a) of Theorem [T] allows us to restrict attention to a strip in which |Im{fc}| is uniformly bounded, 
it then follows that 

- ^ (mW) ^''^ 

as |fc| tends to infinity in this strip. The proof is then complete upon recalling part (b) of Theorem [T] □ 



What do shadows S look like for functions A{-) and S{-) that might arise in applications? Figures [MT^ 
each contain plots of turning point curves T (blue), several of the hyperbola Ec* for various values of C 
(black), the bounds on Im{fc} from from part (a) of Theorem [1] (red), and the shadow region S (green). The 
figures arc all for the potentials (f>{x) = A(x)e*'^'^-'/'^ with A{x) = sech(a;) and S'{x) = sech(x) tanh(x). (This 
example was also considered in the context of the Zakharov-Sliabat spectral problem in [5] and [22].) The 
four figures illustrate the dependence of the shadow region S on the parameter a. (For clarity, these figures 
do not illustrate the e-dependent "collar".) 

Figures [TT] and [T^] each consist of two plots of the same thing; one displayed in the fc-plane and one 
displayed in the A-plane, where A and k arc related by (jl79p as explained in the appendix. When a is small, 
as in these figures, one expects features of the spectrum when viewed in the A-plane to resemble known 
features of the spectrum for the nonself adjoint Zakharov-Shabat spectral problem, which arises by a formal 
limit of a ^ 0. The hyperbolae tend to vertical lines in this limit, and part of the turning point curve 
tends to a fixed curve in the A-plane. Indeed, substituting for k in terms of A from (|179p in the relations 
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Real Pert of k 

Figure 9. The turning point curve 
T and its hyperbolic shadow S for 
a = 0.8944. 
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Real Part of k 

Figure 10. The turning point 
curve T and its hyperbolic shadow 
S for a = 0.6261. 
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Figure 11. The turning point curve T and its hyperbolic shadow S for a = 0.2683. Left: 
the k-plane. Right: the X-plane. 



([53|l defining the turning point curve T, one may select the branch of T for wliich Re{/s} w 1/2 for a small 
and pass to the hmit of a — > with A and x held fixed, to obtain the limiting curve 



lm{A} = ±A{x) , 
Re{A} = -Is'ix). 



(100) 
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Figure 12. The turning point curve T and its hyperbolic shadow S for a — 0.0894. 
the k-plane. Right: the X-plane. 



Left: 



We refer to this limiting curve in the complex A-plane as 7^. Deift, Venakides, and Zhoifl have established 
results analogous to ours in the context of the nonselfadjoint Zakharov-Shabat eigenvalue problem; in these 
results the curve % plays the role of the a-dependent curve T and the family of vertical lines Re{A} = C 
plays the role of the family of hyperbolae Ec . Thus we see that the known "shadow bound" result for the 
Zakharov-Shabat problem can be deduced from our results as a limiting special case. 

The shadow bound for the Zakharov-Shabat eigenvalue problem implies that if 5"'(-) = (so that £"(•) = 
Co for some constant Co), then the eigenvalues lie very close to the real axis or the vertical segment with 
endpoints —Co/2 ± iAmax, ^max := sup^gg A{x), in the complex A-plane. Klaus and Shaw [29] have proved 
that if in addition A(-) is a real function with a single critical point (necessarily a local maximum) then the 
eigenvalues of the nonselfadjoint Zakharov-Shabat spectral problem arc purely imaginary numbers regardless 
of the value of e > 0. Such a result represents a much more precise confinement of the eigenvalues than is 
afforded by the shadow bound with its attendant e-dependent "collar" . It is an open question whether there 
exists such a result in the context of the MNLS spectral problem (|176p . A reasonable guess might be to 
draw an analogy between the vertical "rays" making the shadow in the Zakharov-Shabat problem with the 
hyperbolic "rays" making the shadow in the MNLS spectral problem and thus suppose that if the functions 
A{-) and 5"(-) are such that the turning point curve T is contained within a single hyperbola Ec for some 
C < 0: 

2a^A{xf +aS'{x) ~1 = AC <0, Vx e M (101) 
(this is the MNLS analogue of the condition S'{-) = Co in the Zakharov-Shabat cas(|l, which makes the 
corresponding turning point curve To collapse onto the imaginary axis of the A-planc) , and if an appropriate 
analogue of the Klaus-Shaw critical point condition is satisfied, then the eigenvalues all lie exactly for all 
£ > on the hyperbola E^. Noting that the Klaus-Shaw critical point condition implies that the limiting 
turning point curve 7^ C iM is traced out according to its natural parametrization exactly twice as x increases 
from — oo to +oo, wc might further conjecture that the analogue of this condition for the MNLS problem 
is that T C Ec is covered exactly twice as x varies. But, from ([53|l we see that this occurs for A{-) and 
S"(-) satisfying (jlOip exactly when the amplitude function A{-) is positive and has a single critical point, 
necessarily a local maximum; that is, we expect that the same exact monotonicity condition will be required 
of A{-) as in the Zakharov-Shabat case to obtain exact spectral confinement to a curve, but that the relation 
S'{-)=Co must be replaced by (fTOT|) . 



^This is an unpublished result, see [6], 1321 

^In fact, upon sotting C = {aCo — l)/4 in HlOll l and taking the limit a ■ 
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one obtains exactly the relation S'{-) = Co- 



4. A Family of Special Initial Data 



In this section, we investigate in detail the spectral problem corresponding to the MNLS equation (see 
(|176|) in the appendix) for initial data of the form (j){x^ 0) = A(a;)e*^(^)/'^ where we assume that the amplitude 
and phase functions are given by 



A{x) ~ t^sech(a;) 



and 



S[x) = So + [5 + /itanh(a;)] dx , 



(102) 



where v > 0, 5, i^i, and So are arbitrary real parameters. Following the algorithm described in the appendix, 
our aim is to calculate all of the scattering data for this family. This requires solving the linear system 
(|176p . which we will do by showing that for this special initial data it is essentially a hypergeometric 
equation. Our analysis will be valid for all values of the data parameters v, S, /i, and So and for all 
values of a > 0. More significantly, our analysis will also be valid for all £ > which will allow us to carry 
out a complete semiclassical analysis of the inverse (Riemann-Hilbert) problem in a subsequent publication. 
In connection with this latter application, the family of special initial data under consideration here is 
particularly interesting because it will turn out that for certain choices of parameters in the family and 
certain values of a wc can be in any of the cases where the condition ((2T|) is satisfied for all a; e M, is not 
satisfied for any a; € M, or satisfied for some x and not others. Thus wc have the possibility of providing 
a completely rigorous asymptotic explanation of the effect of the stability threshold on the semiclassical 
dynamics of the MNLS equation. 

The idea of seeking potentials </> for which the spectral problem (|176p might be hypergeometric is an old 
idea for Schrodinger operators. This technique was applied to the nonself adjoint Zakharov-Shabat operator 
(upon which the integrable theory of the focusing NLS equation is based) first by Satsuma and Yajima [37] 
who considered the case (f>{x, 0) = ^'sech(a:), and subsequently by Tovbis and Venakides |38j who generalized 
the technique to potentials of the general forirQ wc consider here. Our work in this section shows that 
potentials of Tovbis- Venakides form are hypergeometric also for the more complicated spectral problem 
(|176p corresponding to the MNLS equation. Moreover, our results generalize those of Tovbis and Venakides 
in the sense that wc recover their formulae by taking a suitable limit corresponding to a ^ 0. 

4.1. Reduction to hypergeometric form. Making the substitution w = e^*'^'^''^''/'^'^ in the linear 
system ()176|) where (f> = A(x)e*'^(^)/^ has the effect of removing the oscillatory terms from the coefficient 
matrix. Thus, with the amplitude and phase functions A and S given by (|102p wc arrive at the following 
equation satisfied by w: 

2il — i/itanh(a;) Aikh'sech{x) 
Aiki^ scch{x) 



2 

dx 



-2n + i/itanh(a:;) 



w. 



(103) 



where 



A 

2 



The key observation we make, in direct analogy with [37j and 
y = tanh(a;) together with the substitution 



a6-l). (104) 
55] , is that the change of independent variable 



a := 



"1 





w = 


'1 








Vi - v'\ 





scch(x) 



w 



(105) 



reduces (|103p to a linear first-order system with rational (in the new independent variable y) coefficients. 
Indeed, a satisfies 



2e(l-2/') 



da 



dy 



2VL 



ifiy 



(106) 



Aikvil-y"^) -217 + (i^ - 2e)?/ 

This 2x2 system has exactly three singular points, y = ±1 and y = oo. Moreover, all three are regular 
singular points, which essentially makes (|106p a hypergeometric differential equation [9]. To make this 
clearer, we follow some steps to convert (jl06p to a more standard form. 



'^Tovbis and Venakides actually assumed that 5 = but elementary symmetries of the Zakharov-Shabat problem corre- 
sponding to the Galilean invariance of the focusing NLS equation imply that their results are easily transformed to completely 
handle the case of 5 7^ as well. 
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The method of Frobenius appHes near each of the (regular singular) points y = ±1. That is, for y in a 
neighborhood of y = ±1 there is a solution a(y) that has a convergent expansion of the form 



n=0 



(107) 



where a„± are vector-valued coefficients. Substituting this series into (|106p and gathering the coefficients 
of like powers of (1 ± y) leads to a hierarchy of algebraic equations relating the coefficients a.n± and the 
exponent p±. The equation arising at leading order is 



T 4ep±ao± = 



2i7 ^ ifi Aikv 

-2n±{i^i-2e) 



ao± • 



(108) 



In order for there to exist nontrivial series solutions (subsequent coefficients vanish for generic exponents 
if ao± — 0) it is therefore necessary that =F4e/5± agree with one of the eigenvalues of the matrix on the 
right-hand side. This is the indicial equation of Frobenius theory. The exponents therefore are 



P± 



ill D, 



(109) 



The case where A:^ G M will be of particular interest shortly, and for such k the two exponents never differ 
by an integer, which means that a full basis of solutions near each of the two singular points y = ±1 may be 
obtained in the form of convergent Frobenius series using both exponents in each case. 

The purpose of calculating the Frobenius exponents is to show that by making the substitution 



(1 + y) 



-(2n-|-iA')/(4e) 



(1-y) 



(2n-,:p)/(4e), 



(110) 



the differential equation p06p is converted into another having the same three singular points, all regular 
singular, but with the property that near each of the points y = ±1 one of the Frobenius exponents has been 
shifted to zero. Thus the equation for b equivalent to that for a via (|110p . 



e{l~y')— = 
dy 





2ikv{\ - y2 



2ikv 
-217 + [i^jL - e)y 



(111) 



is guaranteed to have one solution that is analytic at y = — 1 and another that is analytic near y = +\. The 
zero that has appeared in the coefficient matrix makes it easy to eliminate one of the components and arrive 
at a second-order differential equation; but the process is made even easier by first making the substitution 





(1-'/)- 



leading to the the coupled system 

e— — ~ Zikvc2 
dy 



and 



2iki'ci + [(i/i + e)y - 2Vt]c2 ■ 



(112) 



(113) 



At this point it is trivial to eliminate C2 and obtain a second-order differential equation for ci alone. The 
resulting equation for ci is a standard hypergeometric equation [9] in the independent variable z ~ (1 + j/)/2 
with parameters 



e 

£ 



^ + yi6A^ 



^-/16P^ 



P' 



(114) 



C 



The standard hypergeometric function F{A, B;C; z) is therefore one solution for ci. 
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4.2. Euler transforms and integral representations of Jost solutions. From the point of view of 
the seattcring theory outhned in the appendix, it is necessary to obtain appropriate bases of sohitions of 
the equation (jl76p . equivalent here to the hypergeometric problem (|113|) and use them to compute certain 
Wronskian determinants as functions of k. In principle, all of this information is available in the vast literature 
on special functions (of which and [S] are but two sources). Rather than simply quoting the results we 
need, we proceed directly by deriving integral representations for the required solutions and using these to 
calculate the scattering data. This approach has the advantage that we will have the integral representations 
available for subsequent perturbative analysis of the spectral problem (|176p in the semiclassical limit based 
on Langer transformations. See [7] for an example of this kind of analysis. 

Integral representations for the Jost solutions may be obtained with the help of Euler transforms |25j . 
Suppose ci and C2 have the integral representations: 

c^{y)^ [ C^m-yf'dt and c^iy) = f C^m - vY' dt , (115) 
je Je 

where the contour E, the exponents Pj G C and the functions Cj{ ) (the Euler transforms of Cj(-)) are to be 
chosen so that these expressions satisfy (|113p . Assuming that the contour S is independent of y and noting 
that therefore 

^ = / C^{t){t~yf--Ut, (116) 
dy Jt. 

the first equation of pi3p is satisfied by choosing 

/3i:=/32 + l and C,it) -.^ --^^^C^it) . (117) 

{P2 + 1)£ 

For convenience, we write C{t) = C2{t) and (3 = (32- With these choices, the second equation of (|113[) 
becomes 



+ [{i^i + e)y-2^] / C{t){t~yfdt. (118) 



E 



The aim in satisfying an equation such as this is to equate integrands rather than integrals. As a first step, 
we may of course move the polynomial-in-y coefficients inside of the integrals; in doing this we write y in 
the form y = t — (t — y) so that we may continue to view the integrands as functions of the variables t and 
t — y. Thus we obtain 



e(3{l-t^)C{t){t-yf-^dt+ / [{2ef3 + e + iii)t-2n]C{t){t-yf dt 

E JT. 

C r zLi.2 

£/3 + e + - 



(/3 + l)£ 



C{t){t - yf"-^ dt . (119) 



The next step is to use integration by parts to get all of the exponents of < — ?/ to be the same. Then 
comparison of the integrands will lead to a linear differential equation for C(i), the Euler transform of that 
for C2(y). Since there are three consecutive exponents, two integrations by parts will be required in general 
to achieve this, and we will arrive at a second-order equation for C(i). This would not be advantageous, as 
we began (essentially) with such an equation. We want to reduce the order. 

The point is that by choosing /3 appropriately we may eliminate one of the three consecutive exponents 
in (|119p . a step that will lead to a first-order equation for C(t). Indeed, the quantity in square brackets on 
the right-hand side of (|119|) is a constant depending on /3, and it vanishes if we choose one of the two values 

1 



/3 = -l + — (^-i^i^VAi^^^^^^) ■ (120) 

With either of these two choices for (3 (as long as /3 7^ — 1, and we will discuss this exceptional case once we 
make a proper definition p5ip of /3 below), the right-hand side of pi9p vanishes. Integrating by parts in 
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the first integral on the left-hand side of (|119p . we therefore find that 

[(1 - t^)C{t)\ + [{2el3 + £ + - 2n] C{t) 



{t~yfdt = {), (121) 



assuming the contour E is chosen so that the boundary terms vanish (we will check this later). We solve 
this equation by equating the integrand to zero, yielding the first-order equation for C{t): 

^Jt ~ = + e + i^l)t - 2^] C{t) . (122) 

The Euler transform C{t) is therefore (up to a multiplicative constant) a branch of the multi- valued function 

C{t)^{l^ty+{l + ty'- , where 7± = --^(2e/3 3e ZM T 2f7) . (123) 

2e 

We now suppose that A:^ € M and ji^ — IQk^v^ > 0, which puts k either on the imaginary axis or in 
the real interval (— |/i|/(4i/), |/x|/(4;/)). We arbitrarily make the concrete choice of the + sign in (|120p (and 
interpret the square root there as being positive). The task at hand is to choose the contour S so that 
we obtain integral representations for the Jost solution j^_^\x;k) defined by the chain of transformations 
j^P = vH^wi-^ai-^bi-^c and the boundary condition 



lim jV^^(x-;yfc)e-^^/^ 



(124) 



We therefore attempt to represent the components of j^''(a;; k) = (Jii+(.t; fc), J2i+(a;; k))"^ in the form 

(125) 



J2i+(x;fc) =e-'^(")/"sech(x)e*^°/(2e)gAx/e / c{t){t~yfdt 



We make these formulae concrete in two steps. First we select the branch of C{t) as follows: 

c{t) = Cfi{t-iy+{t + iy- (126) 

where Cq is a constant to be determined, and where the principal branch of the power functions is intended; 
that is, — TT < arg(t ± 1) < tt. Similarly we choose — tt < arg(t ~ y) < tt. Next, we select the contour S 
illustrated in Figure [T^ Our assumptions on k imply that Re{/3} = —1 and Re{7±} = —1/2. The integrals 
in (|125[) are therefore both convergent at t = ±1 for — 1 < y < 1. Moreover, we have 

(1 - t^)Cm -yf = {(i T 1)'/') , as t ^ ±1 , (127) 

so the neglect of boundary terms in the integration by parts argument is indeed justified after the fact in 
this situation. 

To choose the constant Co and verify the boundary condition (|124p we need to calculate the behavior 
of the integrals appearing in (I125P in the limit y ^ 1, corresponding to x +oo. It is easy to see that 
both of these integrals are in fact analytic in a neighborhood of y = 1 , and a simple dominated convergence 
argument shows that 

lim / {t~\y^{t^\y-{t-yf"-^dt^ [ {t - iy++^+^t + 1)'-'- dt , j = o,i. (128) 

(The factor (t — y)^'^-' is uniformly bounded on S independently of y near y = 1, and the remaining factor 
{t — l)''+{t + 1)''^ is in L^{Y,), so indeed dominated convergence applies.) Again recalling Re{7_} = —1/2 
we see that these are both finite under our assumptions on k. This confirms the boundary condition (|124p 
up to the determination of the constant Cq. To get this constant, note that these limiting integrals are 
essentially beta integrals, which is to say that they are ratios of gamma functions. In the case j = 1 the 
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Figure 13. The contour E in the complex t-plane. The branch cuts of the integrand all lie 
on the wavy half-line. 



limiting integrand is integrable at t = 1 as well as i = — 1. Thus the contour E may be collapsed to opposite 
sides of the branch cut giving 

r + l 

{t~iy++^+\t + iy- dt = -2ism{n{j+ + (3 + i)) J {1 - ty++'^+\t + ly- dt 

1 



^_23+''++''-+'5isin(7r(7++/3 + l)) / u'^- {1 - uy++^^+^ du 



= -2 +^++^ ^"i sm(7r(7+ + /? + l))—^, , o , ox 

r(7+ + 7- + P + 3) 



(129) 



Y(-l-7+-/3)r(7++7-+/3 + 3) ' 
where we have used the "beta-gamma" identity 

/\^-i(l-w)^-irf«=M™ , Rc{A}>0, Re{i?}>0 (130) 
and Euler's reflection formula 

sin(7rz)r(z)r(l - z) = tt (131) 
(see [3]). According to this calculation, the boundary condition (|124p requires that we choose 

^ _ (/? + l)£e~^^°/(^^) r(7+ + 7- + /3 + 3)r(-l - 7+ - /3) 
°" 24+7++7-+%fc:/ r(7- + l) ' ^ ' 

and with this choice j+^(a;; k) is determined uniquely. 

This completes our construction of the Jost solution ]''^\x;k). The three remaining Jost solutions may 
also be found explicitly by following similar arguments, but we will be able to calculate all of the scattering 
data from our knowledge of j^'' (x; k) alone. In the case of the specific class of potentials = 4>{x; v, Sq, /i, e) 
under consideration here it is easy to confirm that for each fixed fc G C, whenever v(x; k) satisfies the 
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differential equation (|176p with (j) = 4i{x] v, Sq, fJ-, s), then the vector 



-1 

1 



v(— x; k) 



(133) 



satisfies the same differential equation for the same fc G C, but for a different potential in the same family, 
namely (/)" = (pix; v, —Sq, — /x, e). In this way, we see that, for example. 



i^^\x;k) 



-1 

1 



(134) 



Since our formulae for i^_^\x; k) are valid for all values of the parameters, we also have immediate access to 
integral representations for the components of 'y_ {x;k). 

4.3. Continuous scattering data. We now use the integral representation for ]']^\x]k) to obtain the 
entries of the scattering matrix S(fc) and establish their basic properties. We also will derive formulae for 
the reflection coefficient r(fc) and the reduced reflection coefficient p{z). 

As described in the appendix (see (|184p in particular), the scattering matrix entries S'ii(fc) and S2i{k) 
may be defined in terms of j^i^"* (x; k) by the limits 

5*11 (fc) lim e'^''/^ Jn+{x- k) 



X — > — oo 



and 



52i(fc):= lim e^^''J2i+{x-k). 



X — > — oo 



To calculate S'ii(fc) for fc e M with ^ — 4fc > we must therefore evaluate the limit 



(135) 
(136) 

(137) 



The integral here is not analytic in a neighborhood of ?/ = — 1, but since {t — y)^^^ is uniformly bounded 
for t and y both near —1 for the values of k under consideration, dominated convergence again applies and 
allows us to take the limit under the integral. Thus 



5ii(fc) 



(138) 



Collapsing the contour to the real line and expressing the resulting beta integral in terms of gamma functions 
yields 



^11 (fc) 



r(7_ +(3 + 2) 



(/? + 1)£ r(-7+)r(7+ + 7- + /3 + 3) 

r(7_+/3 + 2)r(-l-7+-/3) 



(139) 



r(-7+)r(7_ + 1) 

Similarly, to calculate S2i{k) for fc^ G M with — Ak^v^ > we must evaluate the limit 



S'2i(/fc) =Coe*^«/(2^) lim 



X — >■ — OO 



*^(^)/" sech(a;)e2^^/'' / {t ~ I)'' + {t + 1)^- {t ~ tanh{x)f dt 



(140) 



Here we cannot immediately pass to the limit in the integral. However, noting that the integrand is integrable 
at < = oo we may first deform the contour S to the top and bottom of the branch cut for real i < — 1, giving 

{t - iy+ {t + I)''- (< - yf dt = 2i sin(7r(7+ + 7_ + /?)) f (1 - t)'^+ (-1 - t)^- {y - tf dt . (141) 

] J —OO 

Next, composing an appropriate y-dependent scaling with a Mobius transformation we make the substitution 



t^-l-il + y) 



1-u 



(142) 
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giving 

f-i 



1 — u H u 



7+ 



i-^-{l-u)-^+-^—'^-'^du. (143) 



Now we may apply a dominated convergence argument to the resulting integral, since Re{7±} — —1/2 and 
Re{/3} = — 1 implies that 



1 ^ 1 + 2/ ■ 
1 — w H u 



and 

Therefore 



lim 



\u^-{l-u)-''+ 



1 ^ 1+y ' 
1 — u ^ u 



-7+-7--/3-2| 



< (1 -u)"°^'^+> = (l-lt)-l/2 



7+ 



Jo 

r(7_ + l)r(-7_-/3-l) 



Since by direct calculation with y = tanh(a;) we have 

e-'^^'^''^" sech(a;)e2^^/^(l + yy+^+'<- 



lim 



we find 



52i(fc) = -Coe-'^°/(2^)2-'3isin(^(7+ +7_ +/?)) 



r(-/3) 

^ g-iSo/e2V/(2e)-n/e+l/2 

r(7_ + i)r(-7„ - /3 - 1) 
r(-/3) 



(144) 
(145) 

(146) 
(147) 



(148) 



= Coe--o/(-')2-^z sin(.(7, + 7- + /? + 3))Ii2L±lE(zI_i^ 

r(-p) 

^ ,ge-.So/s2»M/e r(-7- - /3 - l)r(-7+ - /? - 1) 

2fci. ^''^ ^ r(-7+-7--/3-2)r(-/3) 
^gg-iSo/e2v/£ r(-7_ - /3 - i)r(-7+ - /3 - 1) 
" 2fc^^ r(-7+ - 7- - /? - 2)r(-/3 - 1) ' 

where we used (jl3ip upon substituting for Co from (|132p and the factorial identity r(z + 1) = zr(z). 

All of our analysis so far has assumed that k'^ gM. with — \Qk^v^ > 0. If instead we have k'^ gM. with 
16k^v^ — /i^ > 0, then we claim that all of our results remain valid if we define 



13 



-1 



1 

2e 



2-^2 



IQk^v'^ - > 0. 



(149) 



In this case we have the inequalities Re{/3} < —1 and Re{7±} > —1/2. In fact, all of the preceding arguments 
remain intact exactly as written with the exception of that leading from ()137p through (jl39p and that leading 
from p40p through (|146|) . We emphasize that the results of these calculations remain the same; only their 
proofs must be altered. 

We may combine the formulae for the two sign cases for IGk^v^ — by introducing the function 



R{k) 



-iyjli^ - IQk^v^ , 16A:2j.2 <• ^2 

v/l6fc2i/2 - ^2 ^ 16fc2;y2 > ^2 ^ 

where in each situation the positive square root is meant, and by writing 

/3:=-l- ^(i/i + i?(A:)) , Im{fc2} = 0. 

At this point, we should resolve the question of what happens if 13 ~ —1. From the above definitions it 
follows we may only have /3 = —1 for Im{A:2} = if > 0, and in this case we must have k = 0. When A: = 
the differential equations for ci{y) and C2(y) decouple and become explicitly solvable in closed- form (Eulcr 
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(150) 



(151) 



transforms are not required). One can check that the construction of Jost solutions for = agrees with 
what we have done for fc 7^ in the sense of taking Hmits as fc ^ along the axes. Moreover, the behavior 
at A; = is consistent with what is shown in the appendix (sec p99|) . (|218p . (|223p . and (|224p ) to hold for 
more general potentials. 

The function R{k) has an analytic continuation from the axes Im{fc^} = to the second and fourth 
quadrants, and we define R{k) in the complex fc-planc by introducing appropriate branch cuts as illustrated 
in Figure fT4l Since the quantity R{k)^ = \Qk^v^ — lies in the lower half-plane when Im{fc^} < 0, it can 



Im{/c} 



1 
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R{k) eR+^ R{k) G zM_ 
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Figure 14. The branch cuts of the analytic continuation of R{k) from the axes are shown 
with wavy curves. 



be easily shown by taking the appropriate branch of the square root consistent with the definition of R{k) 
for k in the second quadrant, that for such k we have Re{i?(A:)} > 0. 
In terms of this function, wc then have 

^^2 e 2e )^2 s ~^ 2e ' 

and 

q — -— ^V2^2E^e/-^V2^2e e ^ 

^ V2e 2e /-^ V2e 2e ^ 

The definition of R{k) makes the arguments of all gamma functions in (|152p analytic in k for Imj/c^} < 0. 
This same condition on k then also implies that Rc{l/2 — r2/£±i^/(2e)} = 1/2 — 2Im{fc^}/(£a) > 1/2, which 
verifies the expected fact that 5ii(fc) extends from the axes to an analytic function of k for Im{fc^} < 0. 
While not guaranteed by the general theory, one can nonetheless see directly from (jl52p that 5ii(fc) is 
even in R{k) and hence single-valued and mcromorphic in the whole complex fc-plane (there are poles if 
Im{fc^} > 0). As is consistent with the general theory, S2i{k) docs not have an analytic extension to any of 
the four quadrants; however again we can see from (jl53[) that in this case 6*21 (fc) is even in R{k) and hence 
single- valued in the complex fc-plane. Its only singularities are simple poles lying in all four quadrants. 

33 



Applying the symmetry relations (|230p established in the appendix, we then obtain from (|152p a formula 
for S22{k): 

-I- ii -I- iltWd- -I- iii"! 

^ - r(i + ?-ii(£)i)r(i + ? + ii(£)i) ' ^'''^ 

which extends to an analytic function for Im{/c^} > 0, and from (|153p a formula for 5*12 (/c): 



5l2(fc) 



z£e'^°/-'2-'^/^ r(i - f - g)r(i - i + g) 

-R(fc-')-' N-p/ ijj , 

2e Ml 2e+ 2e ^ (^55) 



^ V 2f: 2.=: /"^ V 2f 2f. > 

Im{fc2} = 0. 



iee*^«/^2-'''/^ r(i - f - f )r(i - i + f ) 



2e 2e 1^ 2e ^ 2e ' 

(The second line follows because for Im(fc^) = 0, either R{k*)* = R{k) or R{k*)* = ~R{k).) The (reflection 
coefficient) function r(fc) = —Si2{k)/S22{k) is therefore 

r{k) = — ^.iAl±^ ^ ' , Im{fc2}=0. (156) 

' ^kiy r(i + f + f) r(-f -^)r(-f + :^) ' ^ ^ ^ ^ 

Since i?(0) = ^^ImIi it is easy to see that r(fc) is an odd function of k that is regular at fc = (in fact in 
this case r{k) is analytic at fc = 0). Therefore the "reduced" reflection coefficient p{z) defined by (|272p is a 
continuous function of z G M given explicitly by 

p{z)^-— 2e ^ /^2+ . 2jn2 + e +2e) I^j^} ^ Q , (157) 

2z. r(i + f + ?) r(-^-£)r(-f + 



where in terms of z = — fc we have 



and 



n = —(4z + l-a6), (158) 

2ia 



Note that for a > the reduced reflection coefficient is never (that is, for no values of the parameters ly, p, 
5, and 5*0 of the potential, nor for any e > 0) identically zero as a function of z S R. This should be strongly 
contrasted with the case a = 0, in which one considers the same initial data in the nonself adjoint Zakharov- 
Shabat problem. For that spectral problem one of the main points made by Satsuma and Yajima [37J is 
that when fj, = there exists a sequence of values of e > tending to zero for which the reflection coefficient 
vanishes identically. This leads to the powerful concept of the "higher-order solitons" — exact solutions with 
simple initial data that can be found by a finite number of algebraic steps as the Riemann-Hilbert problem 
degenerates into a finite-dimensional linear algebra problem. At least in this family of special initial data 
there are no higher-order solitons or reflectionless potentials when a > and therefore the Zakharov-Shabat 
problem must be exchanged for the more complicated problem (|176p . 

4.4. Discrete scattering data. As discussed in the appendix, the discrete spectrum for the problem p76p 
corresponds to the zeros of 5*11 (fc) with Im{fc^} < and those of S'22(fc) with Im{fc^} > 0. By the symmetries 
(|230p and (|23ip it suffices to consider the zeros of S'ii(fc) with Im{fc} > and Rc{fc} < 0. From (|152p one 
can see that zeros of S'ii(fc) occur where the arguments of the gamma functions in the denominator are 
non-positive integers. Since for k in the second quadrant we have both Re{i?(fc)} > and Rc{ri} < 0, only 
one of the gamma functions can contribute zeros of S'ii(fc') in this quadrant. Therefore, the zeros of S'ii(fc) 
with Im{fc} > and Re{fc} < are precisely those values of fc = fc„ such that 

n+lR{k)=(n+l-)e, (160) 



2 ' ' V 2 
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where n is a non- negative integer. It is easy to see that n coincides with 7-|_. Isolating R{k) and squaring, we 
thus find that the desired values k in the second quadrant ar43the solutions with Im{fc} > and Re{fc} < 
of the equation 

+ (iLa + v'^a^)w + — Uv'^a^il - a5) - 4,0^1^ - a^^) , w.^^^ k^ - - [l - a5) , (161) 
16 4 

where L := — e(n+ 1/2) < 0. The discrete spectrum may thus be obtained by solving this quadratic equation 

for w and keeping only those solutions for which Imji/;} < (corresponding to the second quadrant for k). 

This results in a finite number (possibly zero) of roots of Sii{k) in the second quadrant of the complex 

/c-plane. 

If we view L as a general negative real parameter rather than a discrete one, a viewpoint that is increasingly 
accurate in the semiclassical limit £ | 0, we can eliminate the parameter L between the real and imaginary 
parts of (|16ip and hence find a curve C in the (cc, y) = (Re{w}, Im{u>})-plane that is independent of e but that 
necessarily contains the eigenvalues for all e. The eigenvalues for this problem lie exactly on £-independent 
curveQ 

Indeed, in terms of x and y, the real part of (|161|) is 

- y2 _ ^Ly + v''o?x - ^ + \v''o?{\ - Oi8) - ^ = , (162) 

and the imaginary part of (|161[) is linear in L and can be solved to yield 

y Ix^cp-v^ 

L=-- . 163 

a X 

Substitution of (|163p into (|162[) eliminates L and yields the relation 

^2q(^)^!^!^^ q(:^) :=a,2 + ^2^2a;+ l^V(l-a5)~^. (164) 

The curve C defined by (|164p therefore necessarily contains all eigenvalues k in the second quadrant when 
written in terms of w. 

However, not all of the curve C is relevant, because we need to take into account the inequalities y < 
(corresponding to k being in the second quadrant) and L < 0. The latter condition can be translated easily 
into a condition on x with the help of (|163p : since ?/ < and i < we must have x and Ix + oF'v^ being of 
opposite sign, which forces the conditions < x <{). Therefore in addition to lying on the curve C, 

to correspond to an eigenvalue for sufficiently small e, w must also lie in the semi-infinite strip S defined by 

S := {{x, y) G such that y < and -a^v"^ /2 < x < 0} . (165) 

We now claim that the curve C intersects the strip S if and only if ^"^ < 4:^^(1 — aS). Indeed, if this 
condition holds, then (5(0) > 0, in which case the local structure of the curve C near x = y = is a simple 
crossing of two real branches of opposite slopes. Therefore in this situation there is a branch of C near the 
origin with a; < and y < 0, a region that is contained in the strip S. On the other hand, if Q{0) < 0, then 
there are two distinct real branches of C, one with a; > and one with x < —a^i/^, both of which are disjoint 
from S. See Figure [T51 Therefore, if (5(0) < there are no eigenvalues for any e > 0, while if (5(0) > then 
there exist eigenvalues for all sufficiently small e > 0. 

In the context of the particular family of potentials 4> under consideration here, it is interesting and useful 
to compare the threshold condition for creation of eigenvalues with the condition on the parameters v, /i, and 
5 such that the stability condition in pT|) is satisfied for all a; € M. In this specific situation, (|2ip becomes 

a/i tanh(x) + aS + oF'v^ sech^(a-) - 1 > 0. (166) 

Depending on the sign of the minimum of the expression on the left of the inequality in (jl66p occurs either 
as X — > +0O or as a; — > — cx). The minimum value is —a |/Lt| — (1 — a(5). Consequently, if < — (1 — a(5)/a then 



"'^''No superfluous solutions arc introduced because we have already shown that no solutions in the second quadrant can be 
obtained by changing the sign of R{k). 

"'^^Note however, that these curves do not turn out to be hyperbolae of the family Sp. This should perhaps not come as any 
surprise, since for the special initial data under consideration the turning point curve T does not coincide with a hyperbola 
either. See the discussion at the very end of Section |3] 
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Figure 15. Loss of eigenvalues with the variation of S. In all plots, a = 1, /i = 4, and 
1/^2, which makes the threshold for intersection of C and S correspond to 6 = 0. Upper 
left: 5 = —4.1. Upper right: S = —3.9. Lower left: 5 = —2. Lower right: 6—1. 



the condition (|166p is satisfied for all x G M. Furthermore, if > — (1 — aS)/a (in particular ii 1 — aS > 0), 
then the stability condition (|166p fails for some a; e R. 

We may summarize the possibilities for global-in-x stability and existence of discrete spectrum for suffi- 
ciently small e in the following way: 

• 1 — Q!(5 < 0: No eigenvalues exist for any e > 0. 

— If < — (1 — aS)/a then we have modulational stability for all a; G M. 

— If > — (1 — aS)/a modulational instability occurs for some a; € M. 

• 1 — aS > : Regardless of the value of fj., modulational instability occurs for some x e R. 

— If < 2 \iy\ \/l~~aS then eigenvalues exist for sufficiently small e > 0. 

— If > 2 \/l — a6 then eigenvalues do not exist for any e > 0. 

A complete calculation of the discrete scattering data must include the determination of the proportion- 
ality constants 7^ associated with each of the eigenvalues kj in the second quadrant of the complex fc-plane. 
As described in the appendix (see (|227p ). this constant is defined by the relation j^''(a:;; = 7jj^''(a;; fcj) 
holding as an identity for all x G R. Even though it is possible to write down an integral representation 
for j_ (x; kj), it is not convenient to calculate the constants {7^} in this way (with the exception of certain 
problems where it is known that the reflection coefRcicnt vanishes identically; see [57] or [5]). However, 
another way to proceed in special cases such as this is to examine the consequences of the scattering relation 

j^''(x;fc) = 5ii(fc)jL'^(x;fc) + 52i(fc)jL'^(a;;fc), fc^^R. (167) 

It is not generally the case that this equation admits analytic continuation away from the axes, but in the 
special case of the class of potentials under consideration, analytic continuation is indeed possible. By 
continuation along a path beginning on the real fc-axis at a value of k with IQk^pi? > v^, one can show that 
the (nonphysical) continuation of jL^^(.T;fc) into the second quadrant is generically regular at fc = kj. The 
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only poles of this continuation occur if 



'J^ + n = -(m+Ue, (168) 



2 V 2 



where to is a nonnegative integer. This relation also gives poles fc,„ of 6*21 (fc) in the second quadrant lying 
on the hyperbola 

Rc{k„y - Im{/c,„}^ = i(l -a6 + a/i) . (169) 

We call these poles phantom poles. As long as the eigenvalue kj in the second quadrant does not (accidentally) 
coincide with a phantom pole, we learn from (jl67p and the eigenvalue condition Sii{kj ) = that 7^ = 6*21 (kj ) . 
This formula is thus generically valid for the special class of potentials under consideration. 

This expression for the proportionality constants has added significance when interpreted in terms of the 
modified proportionality constants c'^ := Jj / S[i{kj) as one then has 



3 



c" = Res r(fc*)* , (170) 





1 -zp{z) 




1 0" 




1 




p{zr 1. 



zeR. (172) 



where r{k) is the reflection coefficient r(fc) = —Si2{k)/S22{k), or equivalently in terms of the "reduced" 
reflection coefficient, 

-2c," = Jes^p(z*)*. (171) 

Such relations relating the pole data to residues of jump matrix data are indispensible in the analysis of 
Riemann-Hilbert problems of inverse scattering; see, for example [39] . The main idea is to note that the 
jump matrix in the jump condition (j282p for the "reduced" Riemann-Hilbert problem has a factorization: 

"1-Z|p(z)|2 -Zp{zj 

p{zY 1 

The function p(z*)*, with p{z) defined by p57p and mcromorphically continued to the whole complex z- 
plane, has two types of poles in the upper half z-plane: the eigenvalues z„ = — fc^ with A:„ satisfying (|160p . 
and the phantom poles Zm = — with km satisfying (|168p . If the curves supporting the eigenvalue poles 
and the phantom poles are disjoint, then the first and second matrix factors on the right-hand side of 
(|172p may be deformed into the lower and upper half z-plane respectively in such a way as to remove the 
eigenvalue poles from the Riemann-Hilbert problem without introducing any new phantom poles. Unlike in 
the Zakharov-Shabat problem with hypergeometric potential studied by Tovbis and Venakides [551 1 there 
are parameter values in the MNLS case for which the curves supporting the eigenvalue poles and phantom 
poles intersect. Note that if /j, > or /j, < —a^v^jl, the interaction of eigenvalue poles and phantom poles 
is prevented because the hyperbola (jl69p supporting the phantom poles corresponds to a vertical line in the 
w-plane that lies to the right or left, respectively, of the strip S. Phantom poles may only interfere with the 
eigenvalue poles if — q;'^j^^/2 < /i < 0. 

To illustrate these phenomena, we fix the values e ~ 0.075, v ~ 0.6846, 8 ~ 0.5, and \p\ = 0.5 (the discrete 
spectrum only depends on \p\ while the phantom poles are sensitive to the sign of p as well), and allow a 
to vary. According to the criterion established above, eigenvalues will exist for sufficiently small e > only 
for < a < (1 — /i^/(4t/^))/(5 w 1.733. In Figure [TBI we show the emergence of eigenvalues as a is decreased 
below the threshold value of a w 1.733. 

Note that when the hyperbola (|169p that supports the phantom poles passes through the origin in the k- 
plane, that is, when a = l/{6 — p), we also have Q(— (1 — Q!(5)/4) ~ 0, which implies that the curve supporting 
the eigenvalues also passes through the origin in the A:-plane. This situation leads to the interaction of the 
eigenvalues with the phantom poles when p = —0.5, as shown in Figure [T71 

Further decrease of a causes a new phenomenon: the separation of the curve supporting the eigenvalues 
from the imaginary axis. Figure [15] shows the onset of a "spectral gap" between the discrete spectrum 
and the continuous spectrum on the imaginary axis occurring at approximately a = 0.927. The parameter 
a is decreased further in Figure I19[ where we see that the passage through the origin of the hyperbola 
corresponding to the left-hand side of the strip S in the w-plane (this occurs exactly when a = — (i5 ± 
\M^ + 8i^)/(4z^^)) leads to only a small quantitative difference in the eigenvalues and phantom poles and 
the curves that support them. 

Figure [20l shows the final essential transition that occurs as a is decreased to zero: the departure of the 
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Figure 16. The eigenvalues in the complex k -plane are shown with purple dots on the curve 
of the same color that supports them for all e > 0. The green hyperbolae are the images 
of the left and right sides of the strip S in the w-plane. The red hyperbola is the curve 
that supports phantom poles shown as red dots for ^ = —0.5 (that supporting phantom poles 
for jJL = 0.5 is not shown but is necessarily outside the region between the bounding green 
hyperbolae). Left: a = 1.5. Center: a = 1.3. Right: a ~ 1.1. In this figure and all that 
follow in this section, v = 0.6846, b ~ 0.5, = 0.5, and e ~ 0.075. 
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Figure 17. The eigenvalues and phantom poles in the complex k -plane showing onset of 
interaction between eigenvalues and phantom poles. Left: a = 1.004. Center: a = 1. Right: 
a = 0.996. 



hyperbola supporting the phantom poles for /i = —0.5 from the region containing the eigenvalue curve. 
This transition occurs exactly when a = —^/{2v^), which is approximately a = 0.533 for ji = —0.5 and 
v = 0.6846. 

When a becomes very small, one expects some correspondence between our results and those of Tovbis 
and Venakides [38| for the nonselfadjoint Zakharov-Shabat eigenvalue problem which is the formal limit of 
the MNLS spectral problem when viewed in the correct variable A, which is defined in terms of k and a by 
()179l) . Figure [211 shows the discrete spectrum and phantom poles as a tends to zero, both in the complex 
fc-plane and the complex A-plane. In this limit, the discrete spectrum all accumulates near k = ±1/2, and 
the transformation p79p blows up the region near k = 1/2 to reveal a nontrivial limiting structure. This 
limiting structure is in agreement with the results of Tovbis and Venakides [38| in that the limiting discrete 
spectrum consists of equally-spaced points on the vertical line Re{A} = —(5/2 with the phantom poles lying 
on a distinct vertical line: Re{A} = —5/2 + fi/2. 
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Figure 18. The eigenvalues and phantom poles in the complex k-plane showing the appear- 
ance of a spectral gap. Left: a = 0.930. Right: a = 0.924. 
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Figure 19. The eigenvalues and phantom poles in the complex k-plane showing a transition 
in one of the bounding hyperbolae. Left: a = 0.801. Right: a ~ 0.798. 
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Figure 20. The eigenvalues and phantom poles in the complex k-plane showing the phantom 
poles moving away from the eigenvalues as a is decreased. Left: a = 0.7. Center: a = 0.5. 
Right: a ~ 0.3. 
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Figure 21. The eigenvalues and phantom poles in the complex k-plane (top) and X-plane 
(bottom). The red hyperbolae contain the phantom poles for /i = —0.5 (shown with phantom 
pole locations superimposed) and for fi = 0.5 (shown without superimposed phantom poles). 
Left: a = 0.3. Center: a 0.03. Right: a = 0.003. 
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Appendix: Scattering and Inverse Scattering for the MNLS Equation 



Lax pair. Zero-curvature representation. Let e C be a parameter and define 

a \ 4 

Also given a complex-valued function of {x,t) E M?, let 

[A 2ik(j) 
2ikcj)* -A 



and 

r iA2 + 2zP|</.|2 

The simultaneous linear equations 

9v 



-2kk(t> - k£(j)x - 2iak\(j)\'^(f 
-iA2-2ifc2|0|2 



and 



ax 



(173) 



(174) 



(175) 



(176) 



(177) 



for a two-component vector unknown v(a;, i; k) are said to comprise the Lax pair for the MNLS equation ([T]). 
For a general function (j){x,t) the vector v{x,t;k) is overdetermined by the two equations (|176p and (|177p . 
The MNLS equation for (j){x,t) arises as the compatibility condition for the two equations of the Lax pair. 
Indeed, by cross-differentiation one sees that the zero curvature condition 



aL (9B 

e— e- 



[L,B] =0 



(178) 



dt dx 

on the matrices L and B is required for there to exist a full basis of simultaneous solutions of (|176p and 
(|177p . Furthermore, the matrix equation (|178p is equivalent to the MNLS equation governing <j){x,t) as long 
as ^ (the condition (|178p holds trivially if fc = 0). 

While one might think of the MNLS equation ([1]) as a perturbation of the focusing NLS equation when 
a is small, the singular dependence on a in the spectral problem ()176p suggests that at the level of inverse 
scattering we are considering a very singular perturbation. The spectral problem ()176|) can be made to look 
more like the Zakharov-Shabat eigenvalue problem if one introduces the alternate spectral variable A defined 
by the relation 

A=-(2fc-l) or, equivalently, fc = i(aA-fl). (179) 



These relations are clearly singular as a - 
and formally pass to the limit of a 0. 
Zakharov-Shabat problem [IT] : 



> 0. Substituting for k in ()176p we may nonetheless hold A fixed 
The coefficient matrix thus becomes that of the nonselfadjoint 



lim L 

c(->0 
A fixed 



-i\ 



i(f) 
iX 



(180) 



This relation docs not help us to deal with the inverse scattering transform for the MNLS equation, but it 
docs help us to compare some results with known results for the Zakharov-Shabat scattering problem, by 
looking at what happens to the spectrum of p76p for small a in the A-planc rather than the fc-planc. 

Jost solutions and the scattering matrix. Until further notice, we consider t to be fixed and suppress 
all explicit dependence on t. Assume that (pis a complex-valued function of x that satisfies rapidly decreasing 
boundary conditions as x — > ±oo. Then as a; — > ±oo the differential equation (|176p has purely oscillatory 
(rather than exponentially growing and decaying) solutions if and only if k^ € M, and thus the real and 
imaginary axes of the complex fc-plane form the analoguc[3 of continuous spectrum for the differential 
equation (|176[) . 



12 



It is just an analogue because the system is not of the form Cu = ku for any linear operator C; that is, the spectral 



parameter k enters in an essentially nonlinear fashion. 
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(184) 



The Jost solutions of (|176p are, for G R, the columns of two fundamental solution matrices 3±{x;k) 
for this equation having simple asymptotics as x — *■ ±00. By definition for fc^ S M the matrices 3±{x;k) 
satisfy 

e^=:LJ±, lim J±(a;;fc)e-^^'"='/^ =1. (181) 

ox x^±oo 

Since the trace of L is zero, the Wronskian of any two vector solutions of (|176p is independent of x, so by 
taking determinants in the boundary conditions we see that 

det(J±(x;fc)) = 1 (182) 

which shows that if the Jost solution matrices exist for a given k with fe^ € M then they are both fundamental 
solution matrices for (jl76p . Hence they are necessarily related by a constant (independent of x) matrix of 
coefficients: 

3+{x]k) =3_{x;k)S{k). (183) 
The matrix S(fc) for fc^ e M is called the scattering matrix for (jl76p corresponding to the coefficient function 
0. From p82l) and p83)) it follows that 

Sn{k) = det{j^:^\x;k),j^^\x;k)) , 5i2(fc) = det (jf (x; fc), jL'' (x; fc)) 
S2i{k) = det (j^}\x;k),i^^\x;k)) , S22{k) = det (j^}\x; k),i^^\x; k)) 

where (x; k) and j^' {x; k) denote the first and second columns respectively of the matrix J± (x; k). 

Neumann series representations. Analyticity properties. The Jost solutions may be constructed 
from integral equations equivalent to the differential equation and boundary conditions they satisfy. We 
begin by introducing the matrices Y-j-(x; k) related to the Jost matrices J±(a;; k) by 

Y± (x; k) ^ e,-^^'^^'/'^ J± {x; k) . (185) 

By direct calculation using the differential equation (|176p satisfied by 3±{x; k) we find that 

2ik(j){x)e-^^'=/- 
2ik(\)(x)*e^^''l'' 

Also, with the use of the assumption that fc^ e M, the boundary conditions to be satisfied by 3±{x\k) are 
equivalent to the requirement that 

lim Y±(a;;fc) =1. (187) 

X — ^±00 

We therefore seek Y±(x; fc) as solutions of the following integral equations: 

(/)(y)e-2A!// 



dx 



Y± . (186) 



Y±(x;fc)=I+— r 
^ J ±00 



Iterating (|188p once leads to uncoupled scalar integral equations for the four matrix elements of Y-|-(a;; A;). 
In particular, the diagonal elements satisfy 



where the kernel is 



Yii±{x]k)^l + I K{x, z; k)Yii±{z; k) dz 

J±oo 

Y22±{x;k) = 1+ I K{x,z-k*)*Y22±{z;k)dz, 
J ±00 

K{x,z-k):=-—c^{z)* / , 

^ J Z 



(189) 



(t){y)e-'^^^y-'^l' dy. (190) 



Note that the kernel is an even function of the complex variable k. Assuming that (j) G L'^(IR): 

/+00 
\cj){x)\dx<^, (191) 
-00 

we easily get the estimate 

|K(x,z;fc)|<^||0||i.|0(z)|, if sgn(x-^).Im{fc2} >o. (192) 



Seeking to approximate the diagonal matrix elements by iteration starting from the initial guess Yjj± {x; k) = 
1 leads to solutions of (|189p in the form of Neumann series with mth term given by an m-fold integral: 



N,n±{x; k) 



±00 



±00 



K{z2^ zi; k) dzi • • • dz^ 



(193) 



±00 



(or, in the case of Y22±{x; k), N„i±{x; k*)*). We also define No±{x; k) = I. It follows from p92p that 



\N,n-{x;k)\ < 



4|fc|2 



\4>{z)\ dz 



\(l>{z7n^l)\ 



\cj){zi)\ dzi ■ ■ -dZr, 



if Im{fc2} > , 



and likewise that 



\N^+{x;k)\ < 



4ifcr 



|(/)(z)| dz 



if Im{fc2} < 0. 



(194) 



(195) 



By comparison with the exponential series it follows that the Neumann series 

00 

Yii±(x;fc)= ^iV„±(x;fc), ±Im{fc2}<0, 



Tn=0 



(196) 



Y22±{x;k)=Y,N^±i^;k*)*, ±Im{fc2}>0, 



all converge uniformly in compact subsets of the indicated regions of the complex fc-plane and furnish there 
the unique solutions of the Volterra integral equations (|189p . Since the kernel is analytic and even in k it 
then follows by uniform convergence that for each fixed a; £ R all four of these functions are even in k and 
continuous in the indicated regions, also being analytic in the interior. Moreover we have the uniform (L°°) 
estimate 

/4|fc|2 



\Yjj±{x; k)\ < exp 



(197) 



that is valid for j = 1,2 whenever fc^ is in the indicated half-plane of convergence. Together with p92p and 
the assumption that (j> E L^{M.) this estimate shows that 

lim Y:„±(.T;fc) = 1 (198) 

holds for each fixed k in the region of convergence. Now let a; G M be fixed and consider how Yjj±{x; k) 
behaves as fc ^ from each sector of convergence in the fc-plane. Clearly the i°°(M) estimate (|197p is 
uniform for small k as is the bound (|192p of the kernel K. Moreover, since (|192p also shows that A' 
as fc — > pointwise in x and z, a dominated convergence argument applied to the right-hand side of p89p 
shows that 



lim Yj^±(x: k) = I 



(199) 



where the limit is taken from any direction within one of the sectors of convergence. Finally, we consider how 
Yjj± {x; k) behaves as fc ^ 00 within each sector of convergence in the fc-plane. To get the desired estimate, 
we have to note that the apparent fc^ growth of the kernel K should be compensated for by the exponential 
behavior of the integral, so we could first integrate by parts in the definition (|190p to get 



K{x, z; k) 



r(fc2 - 1/4) 



-2A{x-z)/e 



dy 



(200) 



The fraction outside the brackets is now uniformly boimded as fc — )■ 00, and assuming that d(f){x)/dx is in 
L^(E) (which implies that 

II0IIOO :=sup|(/.(x)| (201) 



is finite, so e L°°(M) as well), then we may replace the estimate ()192p by 



\Kix,z;k)\<- 



Ifcp 



£ |fc2 



1/4| 



(2||<^|| 
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if sgn(a; - z) ■ Im{fc^} > . 



(202) 



This bound is uniform in k large. Furthermore, it now follows from dominated convergence that 



lim K{x,z;k) = ~i-\(l){z)f' 



(203) 



with the limit being taken in a strict subsector of a sector for which (|202p holds. Assuming therefore that 
<f) E L-^{R) and that </>' G L-'^(R), we have by dominated convergence that 



lim Nm±ix] k) = 



.a 



\<j>{z)\'dz 



±00 



(204) 



with the limit being taken in the appropriate subsector (note that the presumed conditions on (j) also imply 
that <j) £ L^(R) so the right-hand side is finite). It further follows by dominated convergence applied to the 
infinite Neumann series (|196p that 



lim Yii+(a;; k) — exp [ i— / 



lim Yii^(x; k) = exp 

k — *oo 



+ 00 



\c^{z)\^dz] , 



|0(z)|2dz 



+ 00 



\c^{z)?dz] , 

|0(z)|2dz) , 



lim F22+(x; k) = exp 

k — 'oo 

far 

lim Y22~{x] k) = exp i— / 
fe^oo \ e J_, 

with the limit being taken in an arbitrary strict subsector of the sector of existence in each case. 
Similarly, by iterating p88|) once we find that the functions 

fi2± [x; k) -.^ e2A"/^yi2± (x; k) and f2i± [x; k) -e^^A./^'Faii (x; k) 
satisfy the integral equations 

Yi2±{x;k) = f±{x;k) ^ I K{x, z; k)Yi2±{z; k) dz 

J±oo 

Y2i±ix;k)^f±{x;k*)*+ f K{x, z; k*yY2i±{z; k) dz , 

J±oo 



where the modified kernel is 

k{x, z- k) := K{x, z; fc)e-2'^(^-^)/" 

and where 



2ik 

e 



f±ix;k) 

As before, the assumption that (p E L^(M) gives 

4|fc|2 



4fc2 



z)* / ,/)(y)e-2^('^--)/^d2/. 



dy. 



-2Aiy-x)/e 



±00 



\Kix, z;k)\< ^UWi ■ \(biz)\ , if sgn{x - z) ■ Im{fc^} < , 



and we also have 

|/±(x;fc)| < ^||<^||i, if±Im{A;2}>0, 
estimates that are uniformly valid for small k. The further assumption that 0' G L^(]R) gives 

\K{x, z- k)\ < - ^j^2 _y^ i^Moo + Wh) ■ , if sgn(x - z) ■ Im{k'} < , 

and 



\f±{x;k)\ < I JL_(||^||^ + ||0'||^) , if ±Im{k'} > 0, 



(205) 



(206) 
(207) 

(208) 
(209) 

(210) 
(211) 

(212) 
(213) 



estimates that are uniformly valid for large k. Using (|210p and (|21ip we see that the following Neumann 
series representations converge uniformly for k in compact subsets of the indicated sectors: 



n2±(x;fc)= ^iV™±(x;fc), ±Im{fc2}>0, 

oo 

Y2i±{x;k)=J2N^±i^;k*y, ±Im{fc2}<0, 



(214) 



m=0 



where 



Nm±ix;k):= / K{x,Zm;k) k{zjn, Zm-i; k) 



diCXD 



±00 



K{z2, zi] k)f±{zi] k) dzi - ■ ■ dz„ 



(215) 



The functions Nm±{x;k) are all odd functions of k. Thus, Yi2+{x]k) and Y2i-{x]k) are odd and analytic 
for Im{fc^} > and continuous for Im{fc^} > 0. Similarly, Yi2_(x;fc) and Y2i+{x]k) are odd and analytic 
for Im{fc^} < and continuous for Im{fc^} < 0. If g{x; k) denotes any of these four functions in their sectors 
of continuity, wc have the uniform estimate 



|ff(x;fc)| <^|l</>|li.exp 



A\k\ 



(216) 



A dominated convergence argument applied to (|207p with the help of this estimate then shows that for all 
k in the sectors of existence, 



and similarly that 



lim Yi2±(x;fc) = and lim Y2i±{x;k) = 0, 

X — >±oo X — ^±00 



lim Yi2±{x] fc) = and lim 121^(0:; fc) = 



(217) 



(218) 



with the limits being taken from within the closed sectors of continuity. Also, using (|212p and (|213p we see 
that kg(x] k) is bounded as fc — > cxd within the closed sector of continuity, and by dominated convergence. 



a 1 

lim kNjn±{x]k) = -4>{x) — 
k—>oo I m\ 



a 
i — 

£ 



\4>{z)?dz 



±00 



(219) 



with the limit being taken in any strictly smaller subscctor of ±Im{fc^} > 0. It follows that 

f +00 



lim fcYi2+(x; fc) 



lim fcYi2-_(x; fc) 

k — ^00 



lim kY2i+{x] fc) 

k — >oo 



lim kY2i- (x; fc) 

k — »oo 



\<t){z)\^dz 



X 



2 

a 
2' 

a 



[X ) exp I 



.a 



X 



[X) exp I 



[x] exp 



\(l){z)Ydz 

00 
+ 00 

2 



X 



(220) 



\<P{z)Vdz] , 

|0(z)|2dz 



again with the limit being taken from within a strict subsector of the sector of convergence. 
From this analysis we see that when (jj and </>' are in ij^(R), the Jost solutions 



Yii+{x] fc) 

-Y2l+{X] fc) 



s^^/^ and £\x;k) = 



fi2_(x;fc) 
Y22-{x;k) 



^-Ax/e 



(221) 



-Ax/e 



are analytic for Im{fc^} < and continuous for Im{fc^} < 0. For fixed fc with Im{fc^} < they decay 

exponentially to zero as x — *■ +00 and a; — > — 00 respectively. Similarly, the Jost solutions 

yn_(x;fc) " 
-Y2i-{x;k)_ 

are analytic for Im{fc^} > and continuous for Im{fc^} > 0. For fixed fc with Im{fc^} > they decay 
exponentially to zero as a; — > —00 and x +00 respectively. From (|184|) it then follows that S'ii(fc) 
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j«(x;fc) 



2^^/-' and ]^+\x-k) 



Yi2+(x; fc) 

Y22+{X] fc) 



(222) 



is analytic for Im{fc^} < and continuous for Im{fc^} < while 6*22 (fc) is analytic for Im{fc^} > and 
continuous for Im{fc^} > 0. From (|184p we learn in addition that 



lim Sii(k) = 1 

fc^O 



and 



lim 5*22 (^O = 1 

fc^O 



with the limit being taken from within a closed sector of analyticity in each case. Also 

and lim ^21 (k) = . 



lim 5*12 (fc) = 

fe^O 
Im{fe2}=0 



lim 521 (fc) 

Iin{fc2}=0 



Furthermore, all elements of S(fc) are bounded for k large, and we also have 



lim 5ii(fc) = exp — \ 

k — >oo \ £ 



where 



2 
2 



and 



lim 522 (fc) = cxp 



+00 



|(/)(a;)p dx < 00 



(223) 



(224) 



(225) 



(226) 



and the limit is taken from within a strict subsector of the sector of analyticity in each case. 



Let kj be a zero of 5ii(fc) with Im{fc^} < 0. Since 5ii(fc) is a Wronskian, it follows that the Jost solutions 



j^^\x] kj) and j!:^''(a;; kj) are proportional, with a constant of proportionality that we denote by 7^-: 



:(2) 



3+Hx;kj) = 7jj'L^(.T;fcj) , 



(227) 



Moreover, if we have the strict inequality Im{fc|} < then the left and right hand sides of (|227p both 
represent solutions that decay exponentially to zero as x — > ±00. In other words, these values of k are 
precisely the analogue^ of discrete spectrum for this problem. We further assume that (/) is such that: 

• 5ii(fc) has a finite number of zeros. 

• 5ii(fc) does not vanish for Imjfc^} — 0. 

• All zeros of 5ii(fc) are simple; that is, 5ii(fc) = implies that S'ii{k) 7^ 0. 

Directly from the Neumann series representations, we observe the following symmetries of the Jost solu- 
tions: 



i^i\x;k*)* 



1 
-1 



i'iHx-^k), 



an antiholomorphic symmetry, and 



i^±\x;-k) = 





-1 



'\x;k) 



and 



jf{x;^k) = 



-1 
1 



i^i\x;k) 



a holomorphic symmetry. From the definitions (jl84p we then have 

522(fc) = 5ii(r)* and 52i(fc) = -5i2(fc*)* 

and 

5ii(-fc) = 5ii(fc), 522(-fc) = 522(fc), 5i2(-fc) = -5i2(fc), and 52i(-fc) = -52i(fc) . 



(228) 



(229) 



(230) 



(231) 



Consequently the discrete spectrum for this problem consists of quartets of points of the form (fc, — fc, k* , — fc*), 
distinct because (by assumption) Imjfc^} 7^ 0. The proportionality relation (|227p in conjunction with the 
symmetries of the Jost solutions then implies also that 



j+\x;-kj) = -7jjL^^(a-;-%) 
~j;j'^l\x;k*), 



(2), 



;(2) 



(232) 



j^^^(x;-fc;) = 7;j«(x;-fc*). 



'See the previous footnote. 
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Riemann-Hilbert problem of inverse scattering. Consider the matrix M(fc; a;) defined for x G R in 
terms of tlie Jost solutions as follows: 



M.{k;x) := < 



Sii{k)-+ 

g-A./ej(l)(^.^) 



Im{fc2} < 0, 



Im{fc2} > 0. 



(233) 



It follows from (|184p that det(M(fc; a;)) = 1 at every point in its domain of definition (which implicitly 
excludes the discrete spectrum). Let the (fourfold) discrete spectrum be denoted D. The following other 
properties of M(fc; a;) are fundamental: 

Analyticity: M(fc; x) is analytic for Im{fc^} ^ and k ^ D and takes continuous boundary values 
on the axes Im{fc^} = from each of the four sectors of its analyticity. Moreover, M(A:; x) is 
uniformly bounded for large k. 

Jump Condition: Let M± {k; x) denote the boundary value taken from the region where ±Im{fc^} < 
0: 

u2i _ n (234) 



M±{k;x)-~ lim M{z;x), lm{k^} ^ . 

z — >k 
±Im{z2}<0 

Then the boundary values are related by the formula 

M+{k;x)=M_{k;x)V{k;x), lm{k^} = , 

where 

V{k;x) = e^^'"^/'^Vo(fc)e-^^'^^/'^ , 

and 

^12(fc)^2l(fc) Si2iky 



Vo(fc) 



1 - 



Sll{k)S22{k) S22{k) 

S2i{k) 



1 



(235) 
(236) 

(237) 



Sii{k) 

Note that upon taking into account the symmetries (|230p and (|23ip we sec that in terms of the 
function 

Si2{k) 



r{k) :-- 



S22{k) 



(238) 



±fc2 > 0. 



(239) 



we may write Vo(fc) in the form 

ri±Kfc)P r{k) 

Singularities: The matrix M(fc; x) has simple poles at the points of the finite set D. If kj G D 
with Imjfcj} > and Rc{kj} < 0, then 

0" 



Re^s M{k;x) = lim M{k;x) 



where 



Similarly, 



Cj{x) 



l3 



S'nik,) 



c,{x) 

A, :=A 



^Res M(A:; a;) = ^ linr^ M(k;x 



Res M(fc; a;) = Jim M(fc; x) 



and 



Res^^ M(fc; a;) = ^ lim^^ M(fc; a;) 





Cj{x) 

~c,{xr' 


-c,{xr 





(240) 
(241) 

(242) 
(243) 

(244) 
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Normalization: The matrix M(fc;a;) is normalized in the sense that 

lini M(fc;a;) =1, 



(245) 



with the hmit being taken in any direction. 

Reversing the point of view by taking the function r{k), the fourfold symmetric set D, and the numbers 
c° (which taken together form the set of scattering data associated with (j)) to be given, the above properties 
are said to constitute a Riemann-Hilbert problem for the unknown matrix M(fc;a;). If the matrix M(fc;a;) 
can be determined from the scattering data, then the function (f) may also be recovered therefrom, since from 
the asymptotics in (|205p . ()220p . and (|225p applied to the definition (|233p we see that 



M{k;x) 



2k 



l + o(l) 
</>(a;)* + o{k~ 



a 

2k' 



'{x) + o{k 

1 + 0(1) 

exp 



\<P{z)?dz 







exp 



a 
i — 

e 



\(j){z)\^ dz 



(246) 



as A: — *■ oo in any non-horizontal and non-vertical direction. Consequently, we have the reconstruction 
formula 

2k Mi2{k;x) 



lim 



(247) 



fc^oo a M22(fc; x) 

Note that this Riemann-Hilbert problem differs from others described in the literature even for this specific 
problem (see, for example, [17| ) in that the point of normalization is fc = 0, while the potential cj) is extracted 
via asymptotics of the solution for k large. Indeed, nowhere in the theory of matrix Riemann-Hilbert 
problems is it essential that the point of normalization be any specific value. On the other hand, the fact 
that the potential is extracted by asymptotics at fc = oo is important because while <f> may also be obtained 
from the expansion near fc = 0, these latter formulae involve derivatives of the elements of M(fc; x) with 
respect to x, and controlling these derivatives in addition to the matrix elements themselves in asymptotic 
problems adds an unnecded layer of difficulty. 

Time dependence of scattering data. Now we suppose that 4) is evolving in time subject to the MNLS 
equation ((T|), and we seek to determine the way that the scattering data corresponding to depend on t. 
We need to assume here that for all t we have and 4)x in L^(R) so that the notion of scattering data makes 
sense in the way we have discussed above. We will also assume 4)xx in i^(M). 

Since </> satisfies ([1]), the zero-curvature condition (|178p is satisfied making the Lax pair (|176p and (|177p 
compatible. Therefore, there exists a full basis of simultaneous solutions of (|176p and (|177p . As now 
depends on t as well as x, we now recall the explicit time dependence of the Jost matrices by writing 
3±{x,t]k). Although the columns of J±(a;,i;fc) certainly solve (|176p for each t and although (|176p and 
(jl77p are compatible, the columns of J±(a;,t; k) do not satisfy (|177p because they satisfy time-independent 
boundary conditions as x ^ ±cx). In other words, we can obtain simultaneous solutions of (jl76p and p77p 
by multiplying the columns of J±(a;,t;fc) by scalar factors depending only on t and k and then choosing 
these factors so that (|177p holds for the products. To do choose the factors, it is enough to examine the 
limiting form of the matrix B at ±oo (B has a limit as a; — > ±oo because (/> — > and — > in this limit by 
virtue of the assumption that 0, 0^;, and (l)xx are all in 

e^^(3±{x,t;k)e^^''-^l^] 

These equations hold columnwise throughout the closed sectors of convergence in the complex /c-plane for 
the corresponding Neumann series. 

In particular, both columns of (|248p are meaningful for Im{fc^} = 0. In this case we have the relation 
(|183p which we rewrite here in the form 



for all t). In this way it can be shown that 
B ( J±(a;,i;fc)e'^'*'"^/'^) . (248) 



(j+(.T,<;fc)e*^'*'"-^/'^) = (j_(a;,t;fc)e*^'*'^^/^) (e-'^'*'"-^/^S(fc; t)e''^'*'"-^/^ 



(249) 
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where we have allowed time dependence in the scattering matrix S. Differentiating with respect to t and 
applying (|248p gives 

+ (j_ (a;, t; fc)e'^'*"^/-') • (e-'^'*"^/"S(fe; t)e'^'*"'^^^ . (250) 
Using (|248p again and noting that J_(x,t; k) and e'^^*'^3/e ^j-g invertible we see that 

S(fc; t) = e'^'*'"^/'^S(fc; Oje"*^'*''^/'^ . (251) 
From this relation wc learn the following facts: 

• Sii{k;t) and S22{k;t) are independent of t. Therefore so are the locations of the poles of M(fc;x,i) 
in the complex /c-plane, as well as the derivatives S[i{kj;t). 

• The function r(fc; t) satisfies 

r(fc;t) =r(fc;0)e2'A'*/^ (252) 
The proportionality constants 7^ will also depend on t, and to deduce the time dependence wc first rewrite 
the defining relation (|227p in the form 

i^^^^^^(x,t;k,)) = (e-'^?*/^jL^)(x,t;fc,)) ■ (253) 

Differentiating with respect to t and applying (|248p gives 

B(e'^^*/^j«(x,i;%)) =B(e-'^?*/^jL^)(x,t;fc,)) (^"^^^'/^7,) 

+ (e-^^*/^'jL')(x,t;A:,)) • e| (e^^^^'/^T,) ■ (254) 
Again using (|248p and noting that j^''(x, i; fcj) is not the zero vector, we learn that 

7,(t)=e"2»^^*/^7.(0). (255) 

Given this time evolution of the scattering data wc are led to the following algorithm for solving the 
Cauchy problem for the MNLS equation ([T]). From the initial data (j){x, 0) given in a suitable space, calculate 
the scattering matrix S(fc;0). Check that the function S'ii(A:;0) has a finite number of purely simple zeros 
kj with lm{kj} < 0, and calculate the constants c° := ^j{0)/ S[i{kj;0). Set r{k;0) := -Si2{k;0)/S22{k;0) 
for Im{A:^} = 0, and let D denote the set of points {k £ C|S'ii(fc; 0) = or S22{k] 0) = 0}. The scattering 
data are then used to formulate a Riemann-Hilbert problem: find a 2 x 2 matrix M(fc;x,t), x G M, i > 
such that the following properties are satisfied: 

Analyticity: M(fc;a:,t) is analytic for Im{fc^} 7^ and k ^ D and takes continuous boundary 
values on the axes Im{fc^} ~ from each of the four sectors of its analyticity. Moreover, 'M.{k;x,t) 
is uniformly bounded for large k. 

Jump Condition: Letting 'M.±{k; x,t) denote the boundary value taken from the region where 
±Ini{A:^} < as in (|234p . the boundary values are related by the formula 

M^(fc;x,t) =M_(fc;x,t)e(^^+'A^*)-/^^ ['t^I;Or^'' '^'l 

Singularities: The matrix M(fc;x,<:) has simple poles at the points of the finite set D. If kj E D 
with Im{fcj} > and Re{kj} < 0, then 





-(Aa:+iA^t)o-3/e _|_^2 ^ g 



(256) 



^Res M(fc;x,t) = ^lim M(fc;a;,t) 







and 



where 



Res M(k:x,t)^ lim M(k:x,t) 



c,{x,t) 0_ 
-Cj{x,t)* 
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A, := A 



(257) 

(258) 
(259) 



Normalization: The matrix M(fc;a;,i) is normalized in the sense that 

hm M(fc;a;,t) = I, (260) 

with the hmit being taken in any direction. 
At a given a; G M and i > we say that the Cauchy problem for the MNLS equation ^ has a unique 
solution in the sense of inverse scattering if this Riemann-Hilbert problem has a unique solution for which 
the limit 

2k Mi2ik;x,t) 

<j>{x,t) := \mj ——— (261) 

fc^oo a M22{k\ X, t) 

exists at least in the sense of an approach to fc = oo that is nontangential to the axes Im{fc^} = 0. Lee [3T] 
has generalized the analysis of Beals and Coifman [5] to study the scattering and inverse-scattering maps for 
the derivative NLS equation 

' d^ ^ 2'de ^ '"^W 

with initial data V'(CiO) = "00(0- This problem was shown in the introduction to be related to the MNLS 
equation (IT|) by a Galilean boost with velocity c = a^^. Lee obtains global well-posedness of the Cauchy 
problem in the Schwartz space J5^(M) for a dense subset of initial data, but the theory is not entirely complete 
because certain "nongeneric" initial conditions must be excluded. The nongeneric initial conditions include 
those for which there are "spectral singularities"; that is, the spectral functions analogous to 5'ii(fc) and 
S'22(fc) have zeros on the axes Im{fc^} = 0. 

The solution of the Cauchy problem in the sense of inverse scattering would be expected to agree with 
solutions thereof that have been shown to be unique using Fourier-based (nonlinear) methods. For example, 
in [23j Hayashi and Ozawa considered the Cauchy problem for the derivative NLS equation and showed that 
if ipo is a Schwartz class (^(M)) function, and if 

2 



*7£+o7^+*«7^7(l^lV) = (262) 



!lV^o!|^= / \MOrd^<- (263) 

then there exists a unique global^ solution ■(/'(C, t) that is Schwartz class for all r G M and moreover 
■0 : M ^ is of class C°°(R; 5^). It has also been shown that the constraint (|263|) on the L^(R) norm of the 
initial data may be dispensed with at the expense of knowing existence only for a finite time. For example, 
in [36], Pipolo showed that for initial data t/iq of class H^'°°{E) (meaning that C'^ipoiO^ 'C*'V'o(C)j ^''''PoiO^ 
and ^.'^ip'o'iO s-ll in L^(M) for all integer fc > 0) then there is a number T > such that a unique solution 
■4>{^, r) of ([2B2I) exists of class C°°([-T, T] \ {0}; C°°(]R)). There are many other resuhs in the hterature for 



"'^^It is conjectured in 1231 that the inequahty 12631 1 is sharp at least in the sense that there exist initial data in ,y violating this 
bound leading to solutions that blow up in L°°(R) in finite time. This conjecture is based upon scaling arguments that suggest 
that the derivative NLS equation behaves similarly to the critical collapse case for the one-dimensional semilinear Schrodinger 
equation: 

Despite the Hayashi-Ozawa conjecture, the literature is awash with exact (global) solutions of the derivative NLS equation 
12621 1 that have finite L^(IR) norm but that violate 12631 1. For example, the soliton solutions of 12621 1 given by 

^= — , . \ 02 ' ^U-5o + -Rc{z2}r , e:=-Kc{z^}i + —Kc{z^}T + eo, 

a (ze^ + z*e '■j^ a \ a J a 

where z S C, S M, and G are arbitrary parameters, are all Schwartz class functions uniformly bounded in L^(R) norm 
for each t, but the upper bound exceeds, exactly by a factor of two, the restriction of 12631 1: 



/+00 
mi,r;z)\^ 
- 00 



2n 

sup 

zee . 



It follows also that there exist exact multisoliton solutions (that is, solutions for which r{k; 0) = in the Riemann-Hilbert 
problem), necessarily in with arbitrarily large L^(M) norm. To the degree that approximation by multisoliton solutions is 
possible, there is therefore a dense set of initial data for 126211 with arbitrary L^(IR) norm for which there exists a global Schwartz 
class solution. The validity of the Hayashi-Ozawa conjecture would then seem to rest upon the contribution of nonzero r{k;0) 
to the Riemann-Hilbert problem. Note that although Lee's analysis |31) of the same Cauchy problem using inverse-scattering 
techniques avoids entirely any assumptions about the size of the L-^(R)-norm, it excludes nongeneric data leading to spectral 
singularities and hence is not inconsistent with the Hayashi-Ozawa conjecture. 
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the Cauchy problem associated to the derivative NLS equation (|262p (see in addition [TU], dJ, [12], 
[34| . and [35]) aU differing in detail but all granting existence of a unique solution under either a sufRciently 
small L^(R) norm or a finite lifetime condition. 

All of these results apply equally well to the modified NLS equation ([T]) by a Galilean change of coordinates 
followed by a rescaling of x and t by e. It must be said, however, that these results are unsatisfactory from 
the point of view of the semiclassical (e J, 0) limit for the Cauchy problem for ([T]) with (j){x, 0) = A(x)e''^*^^^/'^, 
because when A and S are fixed the L^(]R) norm on the left-hand side of (|263|) becomes, after rescaling of 
X by e, proportional to 1/e which appears to rule out global solutions. On the other hand, upon rescaling 
of t by e, the fixed lifetime of a solution guaranteed by Pipolo's result becomes arbitrarily small in the limit 
e i 0. So it seems that there are no results in the literature guaranteeing the existence of a unique solution 
to the Cauchy problem for the MNLS equation ^ with initial data of the form 0) = A{x)e'^^^^'^^'^ in the 
vicinity of fixed x and t > when e is sufficiently small. 

Therefore, to deal with semiclassical problems we are forced to take the approach of defining the solution 
of the Cauchy problem by the Riemann-Hilbert problem. Generally speaking, establishing the existence 
of solutions for the Riemann-Hilbert problem is itself a difficult problem requiring careful analysis and the 
correct conditions on the function r(k;0). See [S] or [42j for more details about the corresponding problem 
in a simpler situation, and [31] for information relevant to the equivalent derivative NLS problem. However, 
when one is considering semiclassical asymptotics the problem becomes somewhat simpler because there is 
the possibility of the systematic construction of an explicit parametrix by means of which the Riemann- 
Hilbert problem is reduced to a "small-norm" problem that can be solved uniquely by iteration. This is 
the essence of the steepest descent method of Deift and Zhou. See [57] and [35] for examples of this sort 
of analysis for the case of the Riemann-Hilbert problem corresponding to the focusing NLS equation. We 
intend to carry out such a steepest descent analysis for the Riemann-Hilbert problem corresponding to the 
MNLS equation in a subsequent work. We may have confidence in the success of such a program also 
because certain convergence results for the semiclassical MNLS equation (or more correctly the derivative 
NLS equation) have been obtained by Desjardins, Lin, and Tso using "nonintegrable" methods [16p^ . As the 
effective solution size and lifetime are both large in this situation, these results suggest that global solutions 
indeed may typically exist for large data. However, the approach of [16] fails when the solutions of the 
associated modulation equations become nonclassical, and it is exactly at such junctures that the integrable 
theory typically yields its best results, explaining and giving strong approximations for the wild oscillations 
that are introduced by dispersive regularization. 

Even spectral symmetry. It turns out that our Riemann-Hilbert problem for M(fc;a;,t) has "too much 
symmetry". As we will see below, there is a relation between 'M.{k;x,t) and 'M.{—k;x,t) that will, in 
the course of semiclassical analysis by the steepest descent method, lead to Riemann surfaces with too 
many branch points. For example, in Section [2T3l it was shown that the solutions of the spectral problem 
(|176p corresponding to complex exponential plane-wave solutions of the MNLS equation naturally live on a 
hyperelliptic curve that is two-sheeted covering of the complex /c-plane with four distinct square-root type 
branch points. This surface has genus one, and functions on such a Riemann surface are constructed from 
Riemann theta functions of genus one (that is, elliptic functions). One would therefore expect that the 
potential (f) would itself be constructed from elliptic functions, which appears to contradict our starting 
point: that 4> was of the simple form p^l'^e™^ . 

In this situation, judicious use of factorization identities of Riemann theta functions may be required to 
see that the "true genus" is smaller than first thought. To avoid these difficulties it is desirable to formulate 
a new Riemann-Hilbert problem that effectively quotients out the extra symmetry. This approach was also 
used by Deift, Venakides, and Zhou [2] in their analysis of the Korteweg-de Vries equation, and the specific 
reduction we use here can be inferred also from the paper of Kaup and Newell [28j on the inverse-scattering 
transform for the derivative NLS equation. 

It was pointed out earlier that the Jost functions enjoy a holomorphic symmetry; see ([229]). This symmetry, 
taken together with the implied relation ()23ip and used in the definition (|233p shows that for each t the 
solution M(A:; x, t) of the Riemann-Hilbert problem (assuming existence and uniqueness for the moment) has 



'Although as pointed out in Section |2. 4. ll these authors do not treat the modulationally unstable case. 
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the property that 



M(-fc; X, t) = i^^mik; X, 1)^"'^ = 
It follows that the matrix N(z;x,t) defined by 



Mu{k]x,t) -Mi2{k]x,t) 
~M2i{k;x,t) M22{k;x,t) 



Im{fc2} ^0. 



Mii(fc;x,t) kMi2{k;x,t) 
k-^M2i{k;x,t) M22{k;x,t) 



(264) 



(265) 



is well-defined as a function of z = —k^ for Im{z} ^ 0. From the properties defining M(fc; cc, t) we may then 
deduce a Riemann-Hilbert problem solved by the first row of N(z;x,t). 

It is clear that N(z; a;, t) is analytic for Im{z} ^ 0, and that the boundary values taken for z e R are 
continuous at least for z 7^ 0. To analyze the behavior near z = 0, we see immediately from the normalization 
condition on 'M.{k; x, t) that 



lim Nii{z; x,t) ~ lim N22{z; x,t) ~ 1 



and 



limiVi2(z;a;,i) = 0. 



(266) 



The statement that M2i{k;x,t) — > as /c — > alone is not enough to determine the limiting value of 
N2i{k;x,t); however using the integral equations for Y2i±ik;x,t) it is possible to show that 



lim N2i{z;x,t) 

z— 
±lm{z}>0 



2i 

e 



0(y,i)e~'^^-"^/^""^ dy 



(267) 



±00 



This limit cannot be used to define a normalization condition in the usual way because it has a different 
value in different half-planes for z, and more crucially, it depends on the (unknown) potential (j}{x, t). This is 
the reason why wc will only seek to determine the first row of N(z; x, t) as the solution of a Riemann-Hilbert 
problem (a similar situation occurs in [14]). Being able to find only the first row of N(z;x,i) is not an 
obstruction to determining the potential however, since from the reconstruction formula (|26ip together with 
the antiholomorphic symmetry relations (|228p and (|230p we may equally well find <j){x, t) from the limit 



(t){x,t) 



2 kMi2(k;x,t) 
am. 



2 Ni2{z;x,t) 
lim 



a k^oo Mii{k*;x,t)* a z-»oo Nii{z*;x,t)* 



(268) 



Next, we calculate the jump conditions satisfied by N(z;x,t) on the real z-axis. This is straightforward, 
and the consistency of the calculation (two values of k with Im{fc^} — correspond to each z G M) depends on 
the symmetry relations (|231|) which imply that k~^r{k; 0) is actually a function of z = — fc^ that is bounded 
as z ^ 0. Defining the boundary values by 



N±(z; X, t) 

and noting that in terms of z = — fc^ we have 



lim 'N{w; x, t) , 

w — ^z 
±lm{u)}>0 



z e 



A 



2i 
a 



z + 



one finds that 

N+(z;x,<) = N_(z;a;,i)e(^"+*'^'*''"^/-' 
Therefore, defining 



■l±|r(fc;0)p kr{k-Q) 
±fc-V(fc;0)* 1 



z e 



(269) 
(270) 

±z<0. (271) 
(272) 



p{z) fc-V(fc;0), 

(this is indeed well-defined because r{k;Q) is an odd function of k) we resolve the dichotomy of signs in (|27ip 
and find simply 



N+(z;a;,t) = N_(z;a;,t)e 



l-z|p(z)p 
p(z)* 
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-zp{z) 

1 



z e 



(273) 



Finally, we examine the singularities of N(z; x, t). Near a pole k = kj of M(fc; cc, t) with Im{fc^} < we 



know from (|257|) that 



M(fc;x,t) 



(274) 



"ajCj(fc - kj)-^ + 0(1) flj + 0{k - kj) 
hjCj{k - kj)-^ + 0(1) bj + 0(fc - fcj) 

for some scalar fmictions aj — aj{x,t) and bj — bj{x,t) (and Cj = Cj{x,t) is defined by (|259p ). Applying the 
symmetry ([264]) we see that 



M(fc;x,t) 



'-ajCj{k + kj)-'^ + 0(1) -Gj + 0{k + %)' 
bjCj [k + fc/)-i + 0(1) bj + 0(A: + fc/) 

From these and the definition (|265|) it follows easily that 



and 



N(^;x,t) 



N(z;a;,i) 



ajCj(fc - fcj)"^ + 0(1) kjGj + 0{k - k,j) 
k-\cj{k - k,)-' + 0(1) b, + 0{k - k,) 



-ajCj{k + kj)~^ + 0(1) kjttj + 0(fc + kj) 



k~^bjCj{k + kj)-^ + 0(1) bj + 0{k + k,) 
In turn, these relations are simultaneously equivalent to 

2kjajCj{z - Zj)-'^ + 0(1) kjttj + 0{z - z,) 



N(z;x,t) 
Therefore, we have 



'2bjCj{z - Zj)-^ + 0(1) bj + 0{z - Zj) 
ResN(z;x,t) = lim N{z\x,t) 





-2cj{x,t) oj ' 

a relation equivalent to both relations in (|257|) . In a similar way one derives the condition 



Res N(z; a;, i) = lim N{z;x,t) 



-2z*Cj{x,t)* 




(275) 

(276) 
(277) 

(278) 
(279) 

(280) 



a relation equivalent to both relations in (|258p . 

We arc now in a position to formulate a "less symmetric" Riemann-Hilbert problem to determine (f){x,t). 
Given p{z) for z G M and a set of points zi, . . . , zat in the upper half-plane along with a corresponding set of 
nonzero complex numbers cJ, . . . , c^, seek a two-component row vector {z; x, t) = [ni(z; x, t), n2{z; x, t)] 
with the following properties 

Analyticity: The row vector n'^(z; x, t) is an analytic function of z for z e C\(EU{zi, . . . , z^r, zj^, . . . , z 
that takes continuous boundary values on the real line and is uniformly bounded as z — > oo. 
Jump Condition: With the boundary values defined by 

VLj^{z;x,t) := lim VL'"{w;x,t)^ zGR, 



w — 
±lm{™}>0 



the relation 

n^(z;2;,t) = n^(z; x, t)e(^^+'^'*)^^/-' 



l-z|p(z)|2 -zp{z) 
p{zY 1 



(281) 



(282) 



holds for z G M. 

Singularities: The row vector n-^(z; x, t) has simple poles at the points zi, . . . , zjv and their complex 
conjugates. The residues at these points satisfy 



and 



for J = 1,..., TV. 



Res n"^(z; x, t) — lim VL^{z;x,t) 



Res n"^(z; x, t) = lim n"^(z;a;,i) 

z—z"- z — -fZ* 
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-2c.j{x,t) 

—2z*jCj{x,t)* 




(283) 

(284) 



Normalization: The row vector n'^{z;x,t) satisfies 

limn^(z;a;,i) = [1 O] (285) 

with the limit being taken in any direction. 
From the solution of this Riemann-Hilbert problem, the corresponding solution of the Cauchy problem for 
the MNLS equation ([T]) is then given by 

^M='-lim^^ji^. (286) 
a z^oo ni{z*; x,t)* 

We wish to emphasize that this Riemann-Hilbert problem involves discontinuities only along the real z- 
axis and poles in complex-conjugate pairs. With the extra symmetry removed, this problem is therefore 
quite similar to the Riemann-Hilbert problem associated with the nonself adjoint Zakharov-Shabat eigenvalue 
problem and the focusing nonlinear Schrodinger equation. 
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